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SYMBOLS 


z 

Eleccron,  hole  diffusivity 
Electron,  hole  mobility 
Electron,  hole  lifetime 

*11,  p 

Electron,  hole  diffusion  length  !“(Dn,prn .p^^l 

Jn.p 

Electron,  hole  current  density 

Jnx'Jpx 

X component  of  the  oloctron,  hole  current  density 

JHyJpy 

Y component  of  the  electron,  hole  current  density 

■Eo.Co 

Seturatlon  current  in  Che  base/emitter,  base/collector 

nio>  "I 

Effective  intrinsic  carrier  concentration 
Electron,  hole  concentration 

PO.  n0 

Hole,  electron  electric  field  in  thermal  equilibrium 

^nx.Epx 

X component  of  hole,  electron  electric  field 

Eny.Epy 

Y component  of  hole,  electron  electric  field 

Be-  Ev 

Conduction,  valence  band  edge 

VQF 

Separation  of  quasi -Ferial  potential  <-  vp  - vN) 

NDeff<*> 

Effective  electron  doping  concentration 

**Aeff(*l  Effective  hole  doping  concentration 
NA,ND.Nc  Base,  emitter,  collector  doping  concentrations 


Emitter  length 
Absolute  temperature 

rdy  Bose  de ley  time  in  ! 
Cut-off  frequency 


r corresponding  to  Laplace  transformation 
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REVISED  MODELS  FOR  BIPOLAR  INTEGRATED-CIRCUIT  SIMULATION 


By 

Being  Song  Uu 


systematically  derived  from  the  hole  and  electron  current  and 
continuity  equations.  They  can  be  combined  with  models  now  available 
for  the  p/n  junction  space-charge  layer  to  yield  a network 
representation  for  the  entire  transistor  that  describes  dc.  small- 
signal  and  large-signal  (nonlinear)  transient  responses.  The  models 
consist  of  less  than  ten  lumped  circuit  elements,  yet  their  functional 
dependencies  on  physical  makeup  and  terminal  currents  and  voltages 
yield  greater  accuracy  chan  that  provided  by  previously  available 
models.  The  new  models  apply  for  all  doping  profiles  and  injection 

and  transient  emitter  crowding  at  a higher  level  of  accuracy  than  was 
formerly  attainable.  They  includo  a previously  ignored  energy-gap 
shrinkage  arising  at  high  current  densities  and  Injection  levels  from 


. Traditional  high-level  injection  mi 

Assessment  by  experiments  and  by  computer  simulation  shows  that 
the  new  models  have  advantages  over  their  predecessors  in  simulating 
high  speed  or  high  frequency  or  high  power  bipolar  integrated  circuits. 


CHAPTER  ONE 
INTRODUCTION 

In  Che  development  or  analysis  of  Che  contemporary  bipolar 
integrated  circuits,  an  accurate  and  comprehensive  computer-aided- 
design  tool,  which  includes  process,  device  and  circuit  simulators,  is 
essential  to  Informed  design.  The  utility  of  circuit  simulators 
depends  strongly  on  the  adequacy  of  the  circuit  models  used.  Therefore, 
complete  and  accurate  circuit  models  are  required  for  designing 
integrated  circuits. 

Most  existing  models  for  circuit  simulation,  such  as  Che  Gummel- 
Poon  model  [1],  rely  on  Che  quasi-static  approximation.  This 
approximation  misses  the  detail  of  Che  transient  response  associated 
with  the  time  required  for  carriers  to  propagate  in  any  direction 
within  Che  quasi-neutral  layers.  Thus,  in  the  time  domain  (digital 
circuits),  propagation  delay  is  incompletely  described  by  such  models; 
and  in  the  frequency  domain  (analog  circuits) , the  models  fall  to 
describe  adequately  the  frequency  dependencies  of  the  port  Impedances 
and  the  magnitudes  and  phase  angles  of  the  transfer  functions  ( 2] - ( 17] . 
These  discrepancies  are  particularly  significant  in  circuit  simulation 

are  important  for  1)  transistors  with  thick  base  layers  (power 
transistors) , 2)  transistors  with  small  input  resistances  (emitter- 
couple  logic  circuit),  3)  transistors  in  the  common-base  configuration 


and  4)  heterojunction  transistors  because  the  intrinsic  layers  dominate 
the  transistor  performance  in  these  cases. 

The  purpose  of  this  study  is  to  improve  Che  existing  quasi-static 
models,  such  as  the  Gummel-Poon  model  in  SPICE2  [18),  for  transient  and 

of  a systematic  method  for  the  model  derivation  accomplishes  the 

utilises  an  integration  scheme  to  deal  with  the  time*  and  position- 
dependent  coefficients  of  Che  continuity  and  current  equations  and 

the  models  for  low  and  high  injection  will  Chen  extend  the  model  to 
moderate  injection.  A hierarchy  of  the  models  can  result  either  from 

various  ways  of  partitioning  with  respect  Co  positional  variables  in 
the  quasi-neutral  layers.  This  study  focuses  on  quasi-neutral  layers. 
Appending  models  for  the  junction  space-charge  layers  and  the  extrinsic 


In  bipolar  transistor  modeling  for  circuit  simulation,  it  has  been 
traditional  to  employ  Che  quasi-static  approximation,  thus  neglecting 
the  effects  of  carrier  propagation  delay.  The  prominent  exception  to 
this  statement  has  been  in  small-signal  models  for  frequency  response, 


versus  frequency  representation  of  various  two-port  parameters  such  as 
current  gain,  a or  fi.  This  approach,  however,  does  nothing  to  improve 
the  magnitude  versus  frequency  representation.  In  this  study,  all 
models,  large-signal  and  small-signal,  developed  in  each  chapter,  are 
non-quasi-static  and  accurately  include  carrier  propagation  delay. 

It  has  also  been  traditional  in  transistor  modeling  for  circuit 
simulation  to  assume  idealized  doping  profiles.  Most  of  the  modeling 
details  of  transistor  physics  in  this  study  accommodate  arbitrary 
doping  profiles  and  include  various  physical  mechanisms  in  accordance 
with  Che  following  chapter  outline: 

Chapter  Two — arbitrarily  and  heavily  doped  quasi-neutral  base  and 
emitter  layers  operating  in  low  injection  (including  position- 
dependent  diffusivity,  lifetime  and  energy  gap  narrowing  due  to  dopant- 
carrier  interaction) ; 

Chapter  Three — arbitrarily  doped  quasi-neutral  base  layers  operating  at 
all  injection  levels  (Including  Che  Webster  effect  (19]  and  energy  gap 
narrowing  from  carrier-carrier  interaction  (20]-(23]>; 

Chapter  Four — arbitrarily  and  heavily  doped  quasi-neutral  base  layers 
and  lowly  doped  quasi-neutral  collector  layers  operating  at  all 
injection  levels  (including  Che  Webster  effect,  the  Kirk  effect  (24), 
Fermi-Dirac  statistics,  position-dependent  diffusivity,  lifetime  and 
energy  gap  narrowing  due  to  dopant-carrier  as  well  as  carrier-carrier 

Chapter  Five — exponentially  doped  quasi-neutral  base  layers  operating 
at  all  injection  levels  (including  dc,  ac  and  transient  emitter 


crowding  (25) -(37)). 


Chapters  Two  through  Five  and  concludes  this  work. 

The  models  have  nearly  a traditional  equivalent-circuit  form. 

They  consist  of  connections  of  passive  circuit  elements  and  controlled 
sources,  consistent  with  the  exact  model  for  the  continuity  and  current 
■y  Sah  [38).  The  models  are  not  exactly 

id  negative  capacitances  in  Che  network  representation  for 
transistor.  These  additional  elements  serve  to  Include 
non-quasi-seaeic  effects  while  restraining  the  overall  model 
complexity.  We  incorporate  these  additional  circuit  elements  into 
SPZCE2  by  using  "user-defined  controlled  sources,"  which  are  FORTRAN 
subroutines  accessed  by  SLICE  [39).  Because  the  models  aro  compatible 

SPICE2  by  directly  modifying  the  source  codes. 

developed  provide  better  accuracy  than  does  the  conventional  Cummel- 
Poon  model.  We  assess  Che  improvements  by  comparing  against  experiment 
and  computer  simulations  using  SLICE  and  PISCES  [40)  for  dc,  ac  and 

Semiconductor  Parameter  Analyzer  for  dc  and  s-paramecer  and  the  HP 
4194A  Impedance/Gain-Phase  Analyzer  for  ac  simulations. 

The  analytical  solutions  obtained  during  the  model  derivation 
! co  study  how  device  makeup  and  excitation  affect  the 
r performance.  They  also  provide  information  f< 


CHAPTER  TWO 

ONE- DIMENSIONAL  NON-QUASI-STATIC  MODELS  FOR  ARBITRARILY  AND  HEAVILY 
DOPED  QUASI-NEUTRAL  LAYERS  IN  BIPOLAR  TRANSISTORS 


The  most  widely  used  one-dimensional  model  for  circuit  simulation 
of  bipolar- transistor  circuits  is  the  Integral -charge -control  model  by 
Gummel  and  Poon  [l],  which  Kroemer  (41]  extended  to  characterize 
heterojunction  bipolar  transistors  having  an  arbitrary  doping  profile 

quasi-neutral  base  layer.  To  account  for  time-varying  large-signal  and 
small-signal  excitations  [42] -(43],  it  is  traditional  to  use  the 


quasi-static  approximation,  which  assumes  that  carriers  injected  across 
a junction  prooeed  with  infinite  velocity  during  transients  so  that  the 
steady-state  form  of  the  space  distribution  of  excess  carriers  is 
maintained  at  each  instant  of  time  [44]. 

approximation  while  restraining  model  complexity,  yet  still  accommodate 
arbitrary  doping  and  energy-gap  profiles,  in  both  the  emitter  and  the 
base  quasi-neutral  layers.  Thus  the  models  developed  apply  for  silicon 
and  for  compound  semiconductor  heterojunction  bipolar  transistors.  In 
particular,  we  treat  the  case  of  a heavily  doped  base  layer  in  silicon 
transistors,  which  is  relevant  because  the  small  size  of  individual 
devices  in  an  integrated  circuit  requires  net  doping  levels  above 


approximation 


applications  require  accuracy  in  transfer  functions  for  frequencies  a 
decade  or  more  above  Che  unity-gain  frequency  [45]. 

From  the  viewpoint  of  analysis,  Che  problem  of  developing 
small-signal  and  large-signal  non-quasi-static  models  that  accommodate 
arbitrary  doping  and  arbitrary  energy-gap  profiles  is  challenging. 

Even  for  the  quasineutral  base  layer.  Che  partial  differencial 
equations  encountered  have  time  and  space  dependent  coefficients  and 
are  nonlinear  if  moderate-  or  high-level  Injection  occurs.  Moderate- 
level  injection  in  which  p > 1018  cm'3  occurs  in  transistor  base  layers 
subjected  to  high  current  densities  [46] . In  this  chapter,  we 
eliminate  the  time  dependence  of  the  coefficients  and  the  nonlinearity 
by  restricting  attention  to  the  case  of  low-level  injection  and 

on  the  foundation  laid  here,  we  expand  our  analysis  to  include 
moderate  and  high-level  injection  as  well  as  transient  emitter 

In  Section  2.2,  we  sketch  a systematic  method  for  deriving  non- 
quasi-static  models  for  arbitrarily  and  heavily  doped  base  and  emitter 
layers.  Section  2.3  compares  the  ac  and  transient  responses  of  the 
resulting  models  to  those  of  computer  and  laboratory  experiments.  The 


2.2  Systematic  Derivation  of  the  Model., 

.2.1  Overview  of  the  Derivation 

Wo  will  indicate  here  how  we  include  the  non-quasi-stael 
f propagation  delay  directly  from  the  continuity  equations  a 


neutral  emitter  layers;  as  noted,  we  assume  low  Injection.  The 
derivation  follows  the  scheme  illustrated  in  the  flow  charts  of  Figs. 
2.1(a)  and  2.1(b).  As  Indicated  in  Fig.  2.1(a),  the  continuity 
equations  and  the  current  equations  (which  ore  partial  differential 
equations,  denoted  by  p.  d.  e.)  are  reduced  to  ordinary  differential 

quasi-neutral  layers,  these  ordinary  differential  equations  are 


4<JN< 


:))  f(An(x,t)l  - (l/q)dX( 

- qf*n(x)/(An(x, t) -E  (x, t) 


K.t))/dx 

f qD  (x)d  f(An(x,t)  )/dx 


(2.1) 

(2.2) 


For  low  injection,  the  electric  field  acting  on  the  minority 
carriers  is  approximately  time -independent:  Enix.t)  » Eno ( X ) . Thus 

circuit  models  (47). 

equations  are  obtained  by  integration  and  iteration.  Once  the 
analytical  solutions  emerge,  they  form  the  basis  for  the  development  of 

2.1(b).  The  matrix  elements  ore  infinite  series  in  s,  reflecting  the 
distributed  nature  of  the  quasi-neutral  layers  embodied  in  the  partial 


(CONTINUITY  Eqs,  CURRENT  Eqs.l 

1 M-MM) 


LAPLACE  TRANSFORM  JT fp.d.e.  (*.t)J  «xp(-sl)  dll 


[Nonlinear  o.d.e.(x,s)| 
| LINEARIZE  j 


| Linear  o.d.e.(x.s)  wllh  nonconsl.  coe»Tj 


[integration  and  iteration| 

[analytic  solutions! 

(a) 


The  flow  chart  of  systematic  derivation  of  the  heavily  doped 
non-quasi-statlc  (HDNQSJ  equivalent  circuit  model, 

(a)  Analytic  solutions  including  the  effects  of  carrier 
propagation  delay  derived  directly  from  the  continuity 
foXMd^b’  <b)1C1jCult1reallI,,clon  of  tho  two-port  equation 


differencial  equacions  Chat  describe  them.  To  obtain  an 
equivalent -circuit  model,  which  corresponds  to  ordinary  differential 
equations  in  the  Independent  variable  t,  one  must  truncate  the  matrix 
elements  with  respect  to  s and  take  the  inverse  Laplace  transform. 
Various  truncations  ore  possible.  Particular  truncations  produce 
quasi-static  models,  such  as  the  Gummel-Poon  model.  Ocher  truncations 

elements,  such  as  inductors  or  negative  capacitors.  These  truncations 


2.2.2  Derivation  for  the  Base  Laver 

For  high-speed  or  high-frequency  bipolar  junction  transistors,  the 
submicron  base  layer  has  a peak  doping  concentration  greater  than  lo88 
cm*  to  prevent  punchthrough  [A8]-[50],  According  to  Swirhun  et  al. 
[51],  the  electron  lifetime,  electron  mobility  and  energy  gap  in  the  p- 
type  base  layer  significantly  vary  with  NA(x)  for  NA  > 1087  cm"8.  Thus 
the  coefficients  in  (2.1)  and  (2.2)  dopend  on  x,  which  makes  it 
difficult  to  directly  obtain  analytical  solutions. 

carriers,  we  obtain  Che  electric  field,  En  - d(Ec/q)/dx,  acting  on  the 
conduction-band  electrons 


En(x,t)  » En(x)  = Eno(x)  - (ltT/q)d(ln  »Aaf£)/dX  (2,3) 

where  E0  is  the  conduction-band  edge  and 


where  AEg  Is  Che  apparent  energy- gap  narrowing  defined  In  [51]  which 
includes  Che  effect  of  Fermi  statistics. 

Substituting  (2.3)  into  (2.2)  yields 
XIj,,(x,t>)  - qDn(*)[/(4n(x.c)ld(lnNAoff)/dx  + dX{An(x,t)]/dxI.  (2.4) 


X(An(0,t))  - (nio/NAeff<0)]/lexp(qvBEAT  - 1)1  (2.5) 

X(dn(Ub,t>)  - (njo/Nfleff(Wb) ]X(exp(qvB(./kT  - 1)).  (2.4) 


Ue  substitute  (2.5)  into  (2.1)  and  (2.4)  and  integrate  (2.1)  and  (2.4) 
from  che  emitter  edge  (x  - 0)  of  the  quasi-neutral  base  to  an  arbitrary 
position  x in  the  base.  This  gives 


Z(n'(x,c))  - n0(x)|l  + (Vq)J  X(J,(x,.c)|dxVDn(x')no(x-))  (2.7) 

*(j'(x,t>)  - XfJ'(O.c))  + qj  l2(x',s)Dn(x')X(n'(x',t»dx'  (2.8) 


X(n‘(x,t))  - X(dn(x.t))/X(exp(qvBE/kT)  - 1) 

X{j’(x.t))  - X(jN(x,t))/r(exp(qvBE/kT)  - 1) 

A2(x,s)  - (1  + sro(x))/L2(x). 

Equation  (2.7)  is  of  the  form.  A - a + f)(B) , and  (2.8)  is  of  the 
form  B - b + f2(A).  Thus  A - a + 


fllb  + fj[a  + fx(B)  ] ) . 


- no(x)  [1  + -tlJ'CO.tjljfj  Ajj.jCx.sJ/q  + Jj  B21(x,s>| 


- /(exp<qvgcAT)-l)). 


- £ 1 exp (qvgE/kT) - 1 1 


(2.10) 


If  wo  integrate  (2.1)  and  (2.2)  from  x to  UB  instead  of  from  0 to  x,  wo 
can  obtain  fljn(WB,  t)(  by  following  the  same  procedure.  By  defining 
A.fdutO.t))  - - ijU)  - lgg,  A-/(jN(Mh,t))  - i2  - iCB.  vi  - vBE  and  v2 


A12(s)  1 r f Iexp(qVj/kT) 

A22(s)  I I f(exp(qv2/kT) 


(2.11) 


neglecting  tl 


■ (kT/q)  • A/  JW|> 

-c^c-.c- 

■C*JM 

-C-kC*.* 


» 1 dx  d^  d!t  ,Jub  l d] 


The  results,  (2. l2>-<2.20) , of  this  truncation  have  the  following 
properties.  First,  Ajj  - A21,  as  Is  proved  in  Appendix  A,  which 
confirms  Che  prediction  (521-154]  of  reciprocity  for  bipolar 
transistors  during  transients.  Second,  the  base  delay  time  rd  defined 
by  Elmore  [55]  and  Ashar  [56]  is  rg  1 “ Third,  we  note  Kroemer's 

electron  current  equation  with  charge  control  [41].  His  expression  for 


expression  for  Che  forward  and  reverse  transit  times  obey  Che 
relations:  rcf  - ru  + r21  and  rcr  - r22  + r21. 

2.2.3  Corresponding  Equivalent  Circuits  for  the  Base  Laver 

The  network  representation  for  (2.11)-(2.20)  shown  in  Fig.  2.2 
results  from  redrawing  the  network  implied  by  (2.11)  using  Kirchhoff's 
current  law.  The  element  values  in  Fig.  2.2  follow  from  this  procedure 
or,  more  systematically,  from  employing  the  nonlinear  indefinite 
admittance  matrix  (57).  These  element  values  are 

ijj(s)  - lo[AJ3<s)  + A12(s)JZ(exp(qv ,/kT)  - 1) 

” !„[*(  CjjJAl  * sri2)l-((bxp(lvj/lfr>  -II  3 - 1.  2 

(2.21) 

l12(s>  “ V 'a12(s) 1 f (exp(qv1/kT)  - exp(qv?/kT) ) 

“ 10<1  + sr12)‘1  /(exp(qv1AT)  - exp(qv2/kT)l  (2.22) 

By  caking  inverse  Laplace  transformation  of  (2.21)  and  (2.22),  we 
obtain  the  large-signal  circuit  models  shown  in  Figs.  2.3(a)  and 

RJJ  ‘jj  * (1/ejj  > J0V‘  - Vj  1 - 11  2 <221*> 


CJJ  ■ V ’ll  + t12  H«P<<iyw)-ilAj  ■ Rjj-  '21/Cjj  1 

Ri12(t)  + L di12(t)/dt  - v2  - v2 


(2.21b) 

(2.22a) 


,[exp(qv1/kT)  - exp(qv2/kT)) 


*4= 

Hi(s)( 

) ( 

) *2  2 (S) 

"W77 


Fig.  2.S 


network  representation 


(2.22b) 


non- quasi -Static)  model. 

effects  of  the  quasi-neutral  emitter,  a 


3 also  Include  the 


tent  RjjCjj  i: 


Rjj  i 


present  as  a result  of  the  quosi-stacic  approximation.  In  Fig.  2.3(a), 
the  time  constant  L/R  represents  the  Elmore  or  Ashar  delay  of  the  base 
conduction  current,  a delay  not  included  in  the  Cummel-Poon  model  which 
therefore  has  no  inductance.  The  Cj j in  (2.21b)  times  R in  (2.22b) 
will  give  the  transit  time  rjj  + r12.  From  a physical  viewpoint, 

Rjjcjj  ls  chB  clme  tons tan t for  the  excess  majority  hole  distribution 

minority  electron  distribution  to  rea 
neutral  base  layer.  Since  Che  excess 

excess  electron  concentration  at  each  position  and  time  t< 
quasi-neutrality,  RjjCjj  - L/R.  This  is  consistent  with  (! 

(2.22b). 

forward-active  operation  and  focus  on  the  quasi-neutral  p- 
layer . In  Fig.  2.3(a),  Che  voltage  across  C,i  is  defined 


The  large  - signal  heavily  doped  non-quasi-scaclc  (HDNQS) 
circuit  models. 

(a)  The  inductance  model  for  the  quasi-neutral  base  layer  and 
the  quasi-neutral  emitter  layer  (the  portion  in  the  box); 

(b)  An  alternate  form  of  (a)  for  circuit  application. 


I ‘ VC11 


Dividing  both  sides  of  (2.26)  by  R and  caking  Che  Laplace  transform 


RC11  _ 'tf ■ From  (2.21)  and  (2.22),  we  note  chat 


Substituting  (2.28)  into  (2.27)  yields 


By  taking  Che  inverse  Laplace  transformation  of  (2.29),  we  obtain 


Equation  (2.30)  shows  that  the  HDNQS  model  in  Fig.  2.3(a)  is  equal  to 
chat  in  Fig.  2.3(b)  in  the  forward-active  mode.  Ue  can  prove  Che 
equivalence  of  these  two  models  in  the  reverse -active  and  saturation 
modes  by  following  the  same  procedure.  The  model  of  Fig.  2.3(a) 
reveals  the  physical  insight  regarding  delay  explicitly.  It  is  more 
easily  compared  with  Che  Gummel-Poon  model  and  the  Sah  model  (38)  than 
is  the  model  of  Fig.  2.3(b).  However,  Fig.  2.3(b)  provides  advantages 

common-emitter  or  common-collector  configurations. 

The  circuit  elements  of  the  model  developed  above  are  determined 
by  device  makeup  (geometry,  doping  profile,  etc);  their  accuracy 
depends  on  Che  physical  models  used  for 


energy-gap  narrowing,  mobility 


and  lifetime.  Compared  with 
parameters  of  this  model  chi 

use  it  contains 


le  Gummel-Poon  model,  chere  are  fe’ 
e relevant  physics. 


ill<t) , i22(t)  and  ii2(t)  t 


and  (2.22a)  and  obtain  a non-quasi-s 
i voltage  controlled  current  sources 
it  is  easily  Implemented  in  circuit 


For  the  small-signal  case,  by  assuming  Svj  « kT/q,  we  have 
exp(qVj/kT)  - exp(q(Vj  + «v^  )/kIl  « exp(qVjAT)  • (1  + qiv.  AT)  . 

(2.31) 

Substituting  (2.31)  into  (2.21)  and  (2.22),  we  obtain  the  small-signal 
HDNQS  model  shown  in  Fig.  2.4(a),  in  which 

°JJ  - iy  'jj  ♦ fi2)laxP(qVj/hT)-l)/kT  ; rjj  - rn  /c^  j - 1.  2 

(2.32) 

r - (kX/q)/[Ioexp(qV2AT)];  Lss  -rn-r  (2.33) 

gm  - yaxptqVj/W)  - exp(qV2AT)]/(kT/q)  (2.34) 


where  the  subscript  (ss)  indicates  small-signal  excitation.  Again,  an 
alternate  form  of  the  small-signal  model  is  shown  in  Fig.  2.4(b). 

Since  the  small-signal  circuit  model  is  linear,  it  is  straightforward 
to  prove  that  Fig.  2.4(a)  and  Fig.  2.4(b)  are  equal  by  employing 

If  we  approximate  -Ai2(s)  in  (2.22)  by  1 - sr12,  we  obtain  an 
alternate  HDNQS  model,  in  which  the  series  combination  of  inductance 
and  resistance  in  Fig.  2.3(a)  and  Fig.  2.4(a)  is  replaced  by  a negative 


capacitance  in  parallel 


Bm-  Veil 


Fig.  2.4  Tho  small-signal  heavily  doped  non-quasi-static  (HDNQS) 
circuit  model. 

(a)  The  Inductance  model  for  the  quasi-neutral  base  layer  and 
the  quasi-neutral  emitter  layer  (the  portion  in  the  box); 

Note  that  exchanging  node  E and  node  C gives  another 
Inductance  model;  (b)  An  alternate  form  of  (a)  for  circuit 
application. 


BmVcii 


Fig.  2. 4- -continued 


(2.35) 


The  additional  approximation, 

developed  by  Klose  and  Wieder 
capacitance  model,  is  less  acc 


(1  + Sr^)"1  = 1 - sri2  used  for  this 
in  simulating  the  transconductance  <y2i) 
Lents.  Thus  the  non -quasi -static  model 
[ 7 ] , which  is  equivalent  to  the  negative 
iurate  than  the  inductance  model  [11]. 


Defining  x-0  as  the  boundary  between  the  quasi-neutral  emitter  and 
junction  space-charge  layers  and  x-We  as  Che  emitter  contact  gives  the 
boundary  conditions 

■£(p(0.t)l  - I»i0  /NDef£(0)].t(exp<qvBE/kT)  - 1)  (2.36) 

■£|jp<We.t))  « qSpaX[p(We  ,t)l  (2.37) 

Based  on  these  boundary  conditions,  we  obtain  the  injected  minority 
hole  current  in  the  emitter  layer  by  following  the  method  described  in 
Section  2.2.1 

X(ip(0,t)  = 1^(1  + sret.  /(I  + sred)|  ,f(exp(qvBE/kT)  - 1|.  (2.38) 

emitter  layer,  we  include  recombination  in  deriving  (2.38).  Taking  the 
inverse  Laplace  transformation  of  (2.38)  yields  the  large-signal  model 
for  the  emitter  layer  in  Pig.  2.3,  in  which 


,(0,t) 


(2.39) 


" vbe/ro+  1phqs 


Re  “ vBE''|IEo  'XI>(qVBE/kT)1 

and  lp(jQS.  which  represents  the  non-quasi-static  emitter  storage 
current,  satisfies  the  following  equation 


■*.>£» 


(2.40) 


In  Appendix  B,  we  solve  (2.40)  and  get  iPKQS  in  a form  suitable  for 
circuit  simulation.  By  applying  (2.31)  to  (2.38),  we  obtain  the  small- 
signal  model  for  emitter  layer  shown  in  Fig.  2.4,  where 

re  - <hT/q)/[IEo-exp(qVBEAT)l  , - r^/r^  , r^-  r0(,/ce  (2.41) 


In  Appendix  B,  we  list  the  expressions  for  the  emitter  saturation 
*nt  IEo,  the  emitter  delay  time  re)j  and  the  emitter  transit  time 
By  substituting  IEo  and  rec  into  (2.39)  and  (2.40),  we  find  that 


ii-neutral  quasi-equilibrium  (QNQE)  emitter  model  derived 
[58] . The  additional  resistor  Ree|  in  series  wit 


non-quasi-static  storage  current  in  the  emitter  layer. 

Recently.  Seitchilc  et  al.  [10]  presented  a modified  Gummel-Poon 
model  by  utilising  truncations  of  expressions  for  yjj(s) , y21(s)  and 
0(s) , that  neglect  delay  associated  with  hole  storage  in  the  emitter 
layer.  Comparing  the  delay  times  in  the  base  and  emitter  layers, 


however,  shows  char  the  emitter  delay  time  could  be  ten  times  larger 
than  the  base  delay  time  for  polysilicon  contact  emitters.  Thus  our 

indirect  justification  of  this  statement  in  Figs.  2.10(a)  and  2.10(b). 


2.3.1  Review  of  Previous  Models 

Recent  work  has  sought  to  improve  the  accuracy  of  the  one- 
dimensional Gummel-Poon  model  by  including  non-quasi-statlc  effects. 

In  one  such  effort,  Fossum  and  Veeraraghavan  limited  their  attention  to 
space- independent  base  doping  profiles  in  low  injection  [4).  Klose  and 
Wieder  ( 7 J treated  arbitrary  base  doping  profiles,  using  approximations 
that  yielded  an  accuracy,  for  exponential  doping  profiles  in  the  base, 
equivalent  to  that  of  the  negative-capacitance  model  of  Chen  et  al. 
|H1.  chen  Bt  «1-  also  provided  greater  accuracy  in  what  they  called 
an  inductance  model  for  an  exponentially  doped  profile.  From  now  on  we 
will  refer  to  this  as  the  EDPL  model  (exponential  doping  profile 
inductance  model).  Seitchik  at  al.  [ 10]  developed  a modified  Gummel- 
Poon  model  for  exponentially  doped  profiles  by  using  truncations  of 
expressions  for  yn(s),  y21(s)  and  fit. s). 


2.3,2  Larne-Signal  Models 

To  illustrate  the  improvements  of  the  large-signal  HDNQS  (heavily 
doped  non-quasi-statlc)  models  over  the  original  Gummel-Poon  model  [1], 
we  compare  the  transient  collector  current  predicted  by  the  models  with 
the  exact  solution  in  Fig.  2.5.  The  exact  solution  is  the  solution  of 


Fig.  2.5  The  collector  current  response  to  a large-signal  step  input 
voltage  derived  via  the  Gummel-Poon  (GP)  model,  via  the 
exponentially  doping  profile  inductance  (EDPL)  model  and  via 
the  heavily  doped  non-quasi-statlc  (HDNQS)  model  , compared 
with  the  exact  solution.  A Gaussian  profile  is  used  in  this 
simulation. 


(2.11)  without  truncating  A2i(s).  In  comparison,  we  consider  a p- type 
base  layer  having  a Gaussian  profile  in  which  the  peak  doping  exceeds 
1018cm'3  (heavy  doping).  Wo  focus  on  the  forward  active  mode,  assume  a 
nearly  step  input  voltage  and,  for  the  present,  no  base  resistance. 

In  Fig.  2.5,  we  also  show  the  prediction  of  the  EDPL  model,  which 
uses  space- independent  diffusivity.  mobility,  and  energy  gap  to  get 
compact  analytic  solutions.  If  base  resistance  is  negligible,  the 
delay  exhibited  in  Fig.  2.5  is  entirely  contributed  by  the  intrinsic 
carrier  propagation  delay.  For  an  exact  step-function  voltage,  the 
Gummel-Poon  model  would  show  no  delay  because  it  assumes  Infinite 
minority-carrier  velocity.  The  EDPL  model,  with  delay  Included. 

incorporating  a Gaussian  profile,  having  doping  effects  modeled  after 
Swirhun  (51),  and  delay  gives  a response  within  5%  of  the  exact 

For  simplicity,  we  have  normalized  the  collector  current  by 

we  had  included  changes  in  I0  due  to  the  heavy  doping  effects  and  the 
varying  doping  profile,  the  predictions  of  the  EDPL  model  and  the 
Gummel-Poon  model  would  have  shown  even  larger  derivation  from  the 
exact  solution  than  those  illustrated  in  Fig.  2.5.  Figure  2.6 
illustrates  the  Ashar  delay  time  (56)  versus  the  variance  a2  of  the 
Gaussian  profile  for  the  three  different  models  shown  in  Fig.  2.6.  The 
Gummel-Poon  (GP)  model  shows  no  delay;  for  Che  EDPL  model  Che  delay  is 


substantially  underestimated. 
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Usually  the 


«11-CU>  * P ’tt- 


ideal  step  current -source  Input,  the  transient  response  of  the  non- 
quasi-static  model  will  only  slightly  differ  from  that  of  the  quasi- 
static  model.  However,  an  emitter-coupled  logic  (ECL)  circuit  uses  an 
emitter  follower  as  the  input  stage  to  reduce  the  base  storage  time  and 
to  enhance  the  driving  capability  [59],  Thus  the  input  source 


resistance  for  an  ECL  circuit  is  small  (10-100  ohms).  Because  Che  base 
resistance  of  advanced  transistors  is  also  small  (about  10-100  ohms  for 
the  devices  fabricated  by  MOSAIC  III  (60)).  ECL  circuits  are  voltage- 
driven.  Hence  non-quasi-statlc  models  such  as  those  developed  here  are 
needed  in  principle  for  accurate  ECL  transient  simulations.  The  utility 
of  such  models  is  even  more  pronounced  in  power  and  power-electronic 
applications  where  thick  base  layers  are  used. 


2.3.3  Small-Signal  Models 

2. 3. 3.1  h?1  and  y21 

To  assess  the  validity  of  the  small-signal  model,  we  compare 
several  model  predictions  of  the  phase  of  Che  current  gain  (h2i>  and 
the  phase  of  the  Cransconductance  (y2i>  with  laboratory  experiment  in 
Figs.  2.7(a)  and  2.7(b).  The  experiment  makes  use  of  a large  transistor 
for  which  non-quasi-statlc  effects  become  important  at  low  frequencies 
and  multi-dimensional  effects  are  insignificant  in  low  injection. 
Figures  2.8(a)  and  2.8(b)  illustrate  the  structure  and  doping  profile 
of  the  devices  measured.  We  measured  the  transistors  by  standard  s- 
parameter  techniques.  The  DC  operating  point  is  defined  by  Ic  - 100  pA 
and  by  VBC  - -9.3  volts.  Because  the  doping  profile  is  nearly 
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exponential 


and  the  peak  doping  density  is  less  than  1037  ci'^,  the 
EDPL  model  and  the  HDNQS  model  produce  comparable  results  in  this 
simulation,  and  their  improvements  over  the  GP  model  are  significant. 

To  demonstrate  the  superiority  of  the  HDNQS  model  for  homojunccion 
and  heterojunction  bipolar  transistors,  we  take  a uniform  profile  with 
a linearly  varying  gap  as  an  example  (41 1 and  compare  in  Fig.  2.9  the 
phase  of  y2x(o)  predicted  by  various  models  with  the  exact  solution.  As 
shown  in  Fig.  2.9,  the  HDNQS  model  developed  here  agrees  with  the  exact 
solution  for  frequencies  beyond  the  cutoff  frequency,  whereas  the  EDPL 
model  and  the  GP  model  fail  to  follow  the  exact  solution  even  below 
this  frequency,  A multi-section  EDPL  model  can  perform  comparably  to 
the  HDNQS  model,  but  at  the  expense  of  additional  CPU  time. 

2. 3. 3. 2 yn  and  hu 

In  Figs.  2.10(a)  and  2.10(b),  we  compare  against  exact  solutions 
the  magnitude  and  the  phase  of  y^  predicted  by  the  four  different 
small-signal  models.  The  four  different  models  are  shown  in  Fig.  2.11. 
Figure  2.11(a)  shows  the  small-signal  GP  model  (42),  which  is  based  on 
che  quasi-static  approximation.  Figure  2.11(b)  demonstrates  the 
Seltchik's  model,  which  Includes  Che  non-quasi-stacic  effect  In  the 
base  layer.  Figure  2.11(c)  describes  the  exponentially  doping  profile 
inductance  (EDPL)  model,  which  Includes  che  base  non-quasi-stacic 
effect  by  adding  che  Inductance  into  Fig.  2.11(a)  (6).  Figure  2.11(d) 
exhibits  the  HDNQS  model,  which  Includes  non-quasi-scatic  effects  in 
both  emitter  and  base  layers.  For  simplicity,  we  simulate  a uniformly 
doped  (1017  cm-3)  base  layer  with  0.5  pm  thickness  and  a uniformly 
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rlson  of  yu  predicted  by  four  different  models  shown  in 
..11  with  the  exact  solution. 
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doped  (102®  cm'3)  emitter  layer  with  1.5  pm  thickness  and  with  a 
polysilicon  emitter  contact.  The  discrepancies  between  the  model 
predictions  and  the  exact  solution  shown  in  Figs.  2.10(a)  and  2.10(b) 
imply  that  non-quasi-seatic  effects  are  important  in  the  emitter  layer 
as  well  as  in  the  base  layer.  For  most  transistors  operating  in  low 
injection,  the  non-quasi-scatlc  model  in  Fig.  2.11(b)  provides  accurate 
simulations  of  yn(ft»)  up  to  the  cutoff  frequency  fp.  But  for  those 
transistors  in  which  the  emitter  delay  time  rei  is  comparable  to  /J-rd 
or  “ >ec  » "®®d  the  non-quasi-static  model  in  Fig.  2.11(d)  to 
simulate  yu  adequately  for  f = fT.  Figure  2.12  compares  the  magnitude 
and  the  phase  of  hpp  predicted  by  the  GP  model  and  the  HDNQS  model  with 
experiment.  As  shown  in  Figs.  2.12(a)  and  2.12(b),  the  prediction  of 
the  HDNQS  model  agrees  with  the  experimental  data  for  frequency  up  to 


ten  times  the  cutoff  frequency.  In  contrast,  the  GP  model  fails 
predict  the  experimental  data  beyond  cutoff  frequency. 


A systematic  method  was  developed  to  derive  non-quasi-static 
large-signal  and  small-signal  models  for  quasi-neutral  layers  in 
homojunction  and  heterojunction  bipolar  transistors.  In  this  method,  wt 
analytically  solve  the  continuity  equations  and  the  current  equations. 
To  obtain  the  solution,  we 

• use  Laplace  transformation  with  respect  to  time; 

• integrate  with  respect  to  position; 

• apply  boundary  conditions. 


4000 


static  (HDNQS)  model  and  the  Gummel-Poon  (Cp/model  with  the 
experimental  data. 

(a)  The  magnitude;  (b)  The  phase. 


Fig.  2.1 


inf inice  series  in  5.  Truncation  with  respect  to  s at  different 

simulation  of  high-speed  or  high-frequency  integrated  circuits, 
method  to  Include  moderate-  and  high-level  injection  as  well  as 


CHAPTER  THREE 

ONE-DIMENSIONAL  ALL- INJECTION  NON-QUASI-STATIC  MODELS  FOR  ARBITRARILY 
DOPED  QUASI -NEUTRAL  LAYERS  IN  BIPOLAR  JUNCTION  TRANSISTORS  INCLUDING 
PLASMA- INDUCED  ENERGY-GAP  NARROWING 


3.1  Introduction 

The  purpose  of  this  chapter  is  to  develop  non-quasi-static  circuit 
models  for  quasi-neutral  layers  in  bipolar  junction  transistors  that 
are  general  enough  to  include  all  Injection  levels  and  impurity 
profiles.  These  models,  when  coupled  with  existing  models  for  the 

Cummel-Poon  model  is  based  on  the  quasi-static  approximation,  which 

neutral  layers.  This  assumption  reduces  the  accuracy  of  transient-  and 
frequency-response  simulations. 

Recently  several  non-quasi-static  models  [2] -[12]  that  avoid  the 
inaccuracies  of  Che  quasi-static  approximation  have  been  presented. 

of  applicability  because  either  a specific  baso  doping  profile  or  the 

neutral  layers  under  low  Injection  and  high  injection  based  on  the 
Shockley  equations  and  the  systematic  method  discussed  in  Chapter  Two 
and  [14]  and  reviewed  briefly  here,  a method  chat  applies  for  all 


models  tor  all  levels  of  Injection,  we  use  a fitting  scheme  related  to 
that  In  [61]  to  Include  all  levels  of  Injection  within  quasi-neutral 

simulators  to  optimize  the  fitting  parameters. 

In  treating  high  injection  in  the  base,  it  is  common  to  assume 
that  the  majority  current  density  is  negligible,  which  allows  a simple 
calculation  of  the  electric  field.  In  this  chapter,  we  analytically 

level  of  injection  which  includes  the  special  case  of  low  transistor 


We  also  quantitatively  include  the  effects  of  the  additional 
Coulomb  shielding  that  results  when  the  electron  and  hole  concentration 
exceed  abouc  1017  cm'7,  conditions  that  occur  in  digital  bipolar 
integrated  circuits  [66].  The  presence  of  a high  density  electron-hole 


The  models  developed  have  nearly  a traditional  equivalent-circuit 


controlled  sources,  consistent  with  Sah’s  exact  model  (38]  for  the 
Shockley  equations.  The  models  are  not  exactly  traditional  transistor 
equivalent  circuits  because  they  contain  negative  capacitances  or 
inductances  in  Che  network  representation  for  the  Intrinsic  transistor. 
These  elements  serve  to  include  non-quasi-scaclc  effects  while 
restraining  the  overall  model  complexity.  We  incorporate  these  circuit 
elements  into  SPICE2  by  using  "user-defined  controlled  sources,"  which 
are  FORTRAN  subroutines  accessed  by  SLICE  [39]. 


3.2  Derivation  of  the  Models 

We  follow  the  traditional  practice  of  dividing  Che  transistor  lnc< 
junction  space-charge  and  quasi-neutral  layers.  As  the  main  example 
treated,  we  focus  on  a p-cype  base  layer,  although  the  treatment  con  bi 

One -dimensional  carrier  transport  in  the  quasi-neutral  layers  of 
bipolar  junction  transistors  is  described  by  two  current  equations  and 


jN(x,t)  - q«n(x)n(x,t)En(x,t)  + qDn(xKSn(x,c)/Sx)  (3.1) 
jp(x.c)  - qNp(x)p(x,t)E  (x,t)  - qD  (x)(8p(x,t)/3x)  (3.2) 
Sn(x,t)/8t  - (l/q)(8j[|(x,t)/«x)  - n(x,t)/rn(x)  (3.3) 
8p(x,t)/3c  - -(l/q)(8jp(x,t)/3x)  - p(x,t)/r  (x)  (3.4) 


where  E^p  are  the  electric  fields  acting  on  electrons  and  holes 
respectively. 

In  Che  development  that  follows,  we  begin  in  Sec.  3.2.1  by 


treating  fully  the  simplest 


narrowing  in  Che  p-cype  base  are  boch  negligible.  This  creacmenc 


smaller  than  Che  electron 


current,  and  (b)  the  doping  level  and  Che  hole  and  electron  densities 
are  less  chan  approximately  10*  ? cm'3.  In  Sec.  3.2.2,  we  include  the 

then  indicate  extensions  in  Sec.  3.2.3  to  accommodate  the  large  hole 


E (x,t)  - (j_  + qD  3p/3x)/qa  p - En<x,t). 


(3.5) 


Energy-gap  narrowing  for  low  injection  is  treated  in  the  chapter  two 
and  (14] . In  Section  3.2.3,  we  include  plasma- induced  energy-gap 

Substituting  (3.5)  into  (3.1)  yields  the  nondimensional  relation 


conditions, 

dnj-  (1/2)1  (N^J+ 4nJ(exp(vQFJ/VT)-l)]1/2  - Naj)  j-1,2 (3.7) 

obtained  from  the  quaslneutrallcy  condition.  Here  Vqpj  - vpj  - Vfjj  is 
the  difference  between  che  hole  and  electron  quasi-Fermi  potentials  at 
x - xj . In  (3.7).  Anj  denotes  the  excess  electron  density  at  x - Xj 


For  forward- active  operation,  the  electron  density  at  the  base 
edge  of  Che  collector  depletion  layer  and  the  electron  current  density 
are  correlated,  as  first  explained  by  Middlebrook  [ 62] . This 
interdependence  is  important  particularly  for  Che  analysis  of  base 
pushout.  Thus,  for  full  generality,  use  of  (3.6)  and  (3.7)  for 
forward-active  operation  requires  a simple  Iteration  procedure, 
employing  such  expressions  for  this  interdependence  as  those  found  in 
the  work  of  de  Graaff  (63). 


3.2.1  Method  of  Derivation 

3. 2. 1.1  Arbitrary  doping  profiles  and  low  and  high  injection 
Consider  a p-cype  quasi-neutral  intrinsic  base  layer.  We 

approximate  an  arbitrary  practical  doping  profile  UA(x)  plotted 
logarithmically  versus  x by  a piece-wise  linear  curve  (Fig.  3.1).  In 
Fig.  3.1,  subscripts  J and  j+1  denote  the  edges  (Xj  and  xj+j)  of  Che 
subregion  of  cho  base  that  defines  any  segment  of  the  piece-wise  linear 
curve.  In  the  subregion  Axj  - xj+1  - xj , 

na(*0  - «Aj  exp(qEjX/kT)  ; In  UA(x)  - qEjX/kT  + In  (3.8) 

In  (3.8),  Ej  is  Che  built-in  electric  field  in  an  exponentially  doped 
base  subregion  Axj  for  low-level  injection 

Ej  - - (kT/q)ln(Kj)/AXj  (3.9) 


each  AXj  to  be  thin  enough  that  Dn  is 


through 


thickness. 


illustrate  the  method  of  derivation 


in  Che  simplest  way,  ac  this  point,  we  assume  Chat 
(p/n)((‘p//'n)[lN<x'c)/jp(*.e)l  * 1 i"  eqn.  (3.6).  (This  i 
with  Che  approach  taken  by  many  previous  authors  who  assumed  negligible 
majority  hole  current  in  the  p-type  base.)  We  will  relax  this 
assumption  later,  in  Section  3.2.3. 

Substituting  (3.6)  into  (3.3)  yields 


,["•«  (lnp)/3x  + 3n/SxS(lnp)/3x  + 

i Laplace  transform  of  (3.10)  for  li 
and  high-level  Injection  (p  = dp  s 


• n/r  (x). 

(3 

» Na)  yields  fo: 


d fln/dx  - (l/dXj)lnKj (ddn/dx)  - ( 
where  for  low  injection 


r high  injection 


In  (3.11),  s is  the  Laplace  transformation  variable,  which  equals  d 
Equation  (3.11)  is  an  ordinary  differential  equation  with 
coefficients  independent  of  x.  Combining  its  solution  at  ports  j - 
and  j + 1 - 2 with  the  expression  for  the  injected  electron  curren 
(e  obtain  analytic  solutions  in  a two-port  form 


densicy  in  (3.1), 


n^(s)  -Jv. 

n2(s)/7K 


*JlC) 


layers  In  logarithmic  scale. 


Fig.  3.1 


fitting  arbitrarily  doped 


(AqD/Ax) [ Acoch(A)  + (lnK/2)]/7K 


Aj2  - (AqD/Ax) [Acoth(A)  - (lnK/2)]-,/R 
*12  " ' (AqD/Ax)Acach(A) 

A21  “ ' (*l®/®0*«*eh{*) 

A2  - <lnK/2)2+  Ax2(s/D>. 

Excluding  Che  effect  of  series  resistance,  A^2  - A21.  The  expressions 
for  Che  matrix  elements,  Ajj , shown  above  are  based  on  the 
approximation  chat  the  recombination  lifetime  is  small  compared  with 


To  obtain  the  non-quasl-static  circuit  realization  of  (3.12),  we 
first  truncate  the  matrix  elements  Ajj  to  the  first  order  of  s.  This 
approach  improves  Che  truncation  used  in  [65],  in  which  some  of  the 
coefficients  of  s°  and  sl  terms  are  also  truncated.  The  two-port 
equations  of  (3.12)  correspond  to  either  a hybrld-s  network 
representation  or  a n network  representation  [2],  Here,  to  be 

with  the  form  of  conventional  Cummel-Poon  model,  we  adopt 
isentaclon,  which  yields  Che  large-signal  circuit  model  showi 

^2(3)  + A12(s))An,(s)yK  (3.13] 


» (AqD/Axf2>  (srf/(l  + srd)] 
;(s)  - [Aj2(s)  + A21(s)]An2(s)/VE 


(3.14) 


* (AqD/Axfj)  [sr^d  + srd)  ] At>2(S)/K 
i12(s)  - [-A12(a>]  [AnjtsJ-yic  -2An  <s)/VRj  (3. 

= (AqD/Axfj)  (1  + srd)'1  (An^s)  - An2(s)/K]. 

transit  time,  and  rd  is  the  base  delay  time  defined  by  Elmore  (551  a' 
Ashar  [56],  Expressions  for  these  time  constants  and  for  f i , which 

Appendix  C. 

By  caking  Che  inverse  Laplace  transform  of  (3.13),  (3.1A)  and 
(3.15),  we  obtain  the  large-signal  circuit  model  shown  in  Fig.  3.2  (j 


Jj  JJ 


' SS 


'Jlv—, 


J - 1.  2 


j - (AqD/Axf^)  (rf)[An1(t)/v1l,  Ru~  rd  /Cn 
;2  - (AqD/AxfjK)  (rr)[An2(t)/v2l,  Rjj-  'd  /Cj2 


12(C)/dt  - v,  - v 


R - (v1  - v2  )/(AqD/Axft)  (Anj-  An?/K] , L - r,R. 

The  inductor  L and  the  resistor  Rjj  (R22) . which  ore  not  present  ii 
Gummel-Poon  model,  represent  Che  carrier  propagation  delay  of  the  1 
ent  and  of  the  base  storage  current  respectively.  Ii 
show  an  alternate  form  of  the  large-signal  model  ft 


jl2 

'W %- 


Fig.  3.2 


The  large-signal  equivalent  circuit  models  for  subregion 
AXj.  (a)  The  Inductance  model  for  the  quasi-neutral  base 
layer;  (b)  An  alternate  form  of  (a)  for  circuit  application. 


R * VC11 


Fig.  3.: 


equality  of  Figs.  3.2(a) 


circuit  Implementation.  The 
proved  in  Chapter  Two  and  [14]. 

exp(qVj/mkT)  - exp|q(Vj  + <Vj  )/mkT]  m explqVj/mkT) • (1  + q4Vj/mkT) 

(3. 

where  m-1  for  low  injection  and  m-2  for  high  injection. 

Substituting  (3.18)  into  (3.13),  (3.14)  and  (3.15),  we  obtain  the 
small-signal  modol  shown  in  Fig.  3.3(a),  in  which 


cn  - (AqDn/4xfj)  (rf)(q4N1AT) , ru-  r d /cu  (3.19) 
c22  - (AqD^ixfjK)  (rr)(qAN2/kT),  C^-  r d /c^  (3.20) 
r - (kT/q)/(AqDn/Axf1)  (AN^/K) , L.s  - rdr  (3,21) 
Bm  - <q/kT)/(AqDn/4xf1) [fiNj  - ANj/K)  (3.22) 


where  ANj  and  AN2  denote  the  excess  carrier  densities  on  dc  bias  and 
che  subscript  (ss)  indicates  small-signal  excitation.  An  alternate 
form  of  the  small-signal  model  appears  in  Fig.  3.3(b). 

The  models  with  the  inductance  (Figs.  3.2a  and  3.3a)  reveal  delay 
explicitly.  They  are  more  easily  compared  with  the  Gummel-Poon  model 

models  with  the  voltage-controlled  current  sources  (Figs.  3.2b  and 
3.3b),  however,  provide  advantages  in  circuit  implementation. 
Approximating  -A12(s)  in  (3.15)  by  1 - srd  gives  the  negative- 

combination  of  inductance  and  resistance  in  Fig.  3.2  and 
replaced  with  a negative  capacitance  in  parallel  with  a c 


Fig.  3.: 


Fig.  3.3  The  small-signal  equivalent  circuit  model  for  subregion  Ax*. 

(a)  The  inductance  model  for  the  quasi-neutral  base  layer; 
Note  that  exchanging  node  E and  node  C gives  another 
inductance  model;  (b)  An  alternate  form  of  (a)  for  circuit 
application. 


9mvc1l 


Fig.  3. 3- 'continued 


- 1/R. 


- 1/r 


(3.23) 


The  additional  assumption,  (1  + srj)'1  = 1 - srd.  makes  the  negative- 
capacitance  model  less  accurate  than  the  inductance  model  in  simulating 
transconductance  <y2i)  at  high  frequencies  and  in  representing  the 
large-signal  transient  response  (11). 


3- 2. 1.2  Equivalent -Circuit  Models  for  ALL  Injection  Levels 

Since  an  analytic  model  for  the  base  layer  in  moderate  injection 
is  difficult  to  get  because  the  continuity  equation  is  nonlinear,  we 
fit  parameters  (such  as  rt,  rc  and  rd  and  the  forward  and  reverse 
collector  currents.  If  and  Ir)  based  on  the  low-injection  and  high- 
ie  model  for  all  injection  levels  |61], 


s injection  level  gives  the  gene 

r L,  Cji  (C22)  and  Rlt  (R22)  c. 
rsus  injection  levels.  The 


For  example,  fitting  If  and  I, 
expression  for  R in  Fig.  3.2(a)  for  a 
determined,  the  general  expressions  f. 
be  deduced  by  fitting  rd,  rf  and  r ~ vi 
resulting  models  for  all  injection  levels  have  the  same  form  as  Fig. 
3.2  and  Fig.  3.3,  though  the  expressions  for  the  circuit  elements 
differ  from  those  applying  for  low  and  high  injection. 

For  the  all-injection  large-signal  model,  we  denote  the  forward- 
active  transit  time  by  TF,  reserving  tf  for  low  injection  and  high 
injection.  Similarly  TDF  replaces  rdf,  where  the  f refers  to  forward 
active  operation,  IF  is  a fitting  parameter  related  to  the  forward 
current,  etc.  The  circuit  elements  of  this  model  are  defined  by 
Cu  - TF-IF-IAn^O/Vj],  Ru-  TDF/Cn  (3.24) 


C22  - TR  IR  (an2(t)A2l.  R22-  TDR/C22  (3.25) 
R - (Vj  - v2  )/[IP-ani-  IR-irij]  (3.26) 
L " <V1  ' v2  )/I1F-4n!  /TOf  - IR-inj/TDR)  (3.27) 

cll  - IF-If-IqiNj/CT),  ru-  TDF/CU  (3.28) 
c22  - TR  Ir-lqlSNjAT],  r^-  TDR/c22  (3.29) 
r - kT/(q!r-6N2),  Lss  - TDR/r  (3,30) 
6m  - Idf-AN!  - Ir-ANgJ/kT.  (3,31) 


and  TDR  are  listed  in  Appendix  C. 

We  have  assumed  that  the  functional  dependencies  on  current  of  time 
constants  (TF,  TR.  TDF  and  TDR)  of  the  large-signal  model  are  the  same 
as  those  of  the  small-signal  model  at  all  injection  levels.  This 
implies  Qg(base  storage  charge)/Ic  - dQg/dIc  [42) . This  assumption 
introduces  no  error  in  low  injection  and  high  injection.  For  moderate 
Injection,  we  retain  this  assumption  to  minimize  the  number  of  model 
parameters.  To  achieve  better  accuracy,  we  may  follow  the  same  fitting 


scheme  to  fit  four  more  parameters  for  the  small-signal  model. 

Because  of  the  reciprocity  shown  in  (3.12),  the  forward  delay  time 
TDF  is  enough  in  principle  to  characterize  the  delay  in  the  forward- 
active  mode  as  well  as  in  the  reverse-active  mode.  In  any  subregion  Ax 
- (Xj  - x2).  however,  we  may  have  high  Injection  at  xt  and  low 
injection  at  x2.  Then  reciprocity  is  invalid  and  TDF  fails  to  describe 
the  delay.  Thus  we  need  to  fit  an  additional  time  constant  TDR,  which 


Is  controlled  by  the  Injection  level  at  x2.  to  achieve  better  accuracy. 
We  can  also  replace  Anj/NAj  (see  Appendix  C) , which  corresponds  to  the 


+ Axj/2 

I (Anj/NAj)dx/(AXj/2  + Axj  _2/2) 

Axj.i/2 

to  reflect  the  Injection  level  in  subregion  Ax  j . Although  we 
benefitted  from  the  insight  in  [61],  their  fitting  procedure  ignored 
this  possibility.  Therefore,  their  prediction  of  the  reverse  collector 
current  and  the  reverse  transit  time  differs  from  our  prediction. 

The  number  of  the  subregions  chosen  In  Che  models  is  mainly 
determined  by  doping  profile  and  accuracy  required.  In  the  simulation 
of  heavily  doped  layers,  the  position  dependence  of  diffusivlty, 
lifetime  and  energy  gap  also  need  to  be  considered  while  partitioning 
base  layers  because  positional  averages  of  these  parameters  are  used  in 

The  models  developed  here  are  in  the  common-base  configuration.  To 
derive  the  common-emitter  or  the  common-collector  configuration  of 


hoose  the  hybrid-* 
se  the  models  showr 


f exploit  the  Indefinite  admittance  matrix 
lecwork  representation  in  the  realization  c 
direct  way  to  perform  this  transformation 


[57]  o 


S*  5.2(b)  and  Fig.  3.3(1 
show  a two-section  small-signal  model  in  Fig.  3.4 
emitter  configuration.  To  implement  the  models  ii 
simulators,  such  as  SPICE,  we  employ  the  method  o: 
Controlled  Source  [39]  t 


a circuit 
Jser  Defined 


nonlinear  circuit  elements 


n Appendix  C,  we 


show  three  volcage  controlled 
the  large-signal  model. 


3-2.2  Inclusion  Of  Plasma -Induced  Energy-Cap  Narrowing 

The  concept  of  apparent  energy- gap  narrowing  AEg  Is  presently  widely 
used  [66).  (67).  Apparent  energy-gap  narrowing  is  defined  by  the 
relation,  pn  - n1exp(AEg/kT) , which  thus  includes  the  effects  of  Fermi 
statistics  present  in  the  more  physically  correct  relation,  pn  - 
niexP(AEg/kT)  Fl/2(<qv„  - Ec>/kT|exp[-(qvN  . Ee)/U] Fygl  (By  - 
qvP)/kT)exp(-(Ev  - qvpJ/kT)  where  AEg  is  the  true  energy-gap  narrowing 
measured  between  the  conduction-band  edge,  Ec,  and  valence-band  edge, 
Ev,  and  ”hero  f1/2  1s  cha  Fermi  integral  of  order  1/2.  In  this  chapter, 
we  will  follow  the  prevailing  practice  of  using  the  apparent  energy-gap 
narrowing,  in  which  the  parameter  AEg  includes  the  effect  of  Fermi 


Energy-gap  narrowing  6Eg  due  to  high  concentrations  An  of  excess 
carriers  in  lightly  doped  base  layers  was  measured  in  (20).  This 
plasma- induced  band  gap  narrowing  causes  the  electric  field  acting  oi 
holes  to  differ  from  chat  acting  on  electrons  by  (l/q)8AEg/3x  [68]. 
Including  this  term  in  (3.5)  and  deriving  the  electron  current  under 
the  assumption.  <p/n)  ((ip/Pn)  (JN(x,  t)/jp(x,  C) ) . 1.  gives 

- qDn(28An/3x  - (l/kDAnflAB  /Sx|  . (3.3 


Equation  (3.32)  is  a nonlinear  equation  because  AE„  depends  on  An 
determined  in  [20] -[22] 


(plasma)  m a0(kT/q)ln(An/ni)  - aj , 


1.33) 
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Fig,  3.4  The  tvo-section  smell-signal  equivalenc  circuit  model  for  a 
Intrinsic  base  layer  in  the  forward  active  mode. 


which  shows  the  dependence  of  AEg(plasma)  on  Che  chemical  pocential 
<kT/q)ln(dn/n£) . In  (3,33),  nj  is  che  inecinsic  carrier  concencracion. 
Subscicucing  (3,33)  inco  (3.32)  yields 


j„(*.c)  - qDn(2  - a0/q)San/dx  - qD  f 3An/«x  (3.34) 

According  Co  daCa  measured  in  [20],  oq  is  approximaCely  0.7-q  and  aj  = 
265  mev  for  An  5 SxlO17  cm'3,  and  a0  = q and  aj  » 414  mev  for  An  > 
5xl033  cm'3.  Thus,  for  An  £ SxlO33  cm*3,  fp  is  abouC  1.3  inscead  of  2. 
Thac  implies  thac  Che  base  cransic  cime  including  AEg  in  high  injocclon 
will  be  abouC  1.5  cimes  larger  chan  char  excluding  AEg  in  high 
injecclon.  We  can  include  chis  plasma  effect  in  Che  models  developed  by 
adding  fp  Co  Che  ficcing  parameCers  in  Appendix  C. 


We  now  relax  Che  assumpeion,  (p/n) (ap/p„) ( jN(x, C)/jp(x, C) ] . 1,  made 

For  cransiscors  operacing  ac  high  currenc  levels,  che  minoricy 
eleccron  currenc  in  che  p-cype  quasi-neucral  base  layer  increases  less 
scrongly  wich  vBE  chan  does  che  majority  currenc,  because  of  excrinsic 
resiscances  or  base  pushouc.  Thus,  even  for  forward-active  operation, 
che  majoricy  hole  current  could  be  comparable  Co  che  minoricy  eleccron 
currenc  ac  high  injecclon  levels.  Furthermore,  for  cransiscors 
operacing  in  che  reverse-active  or  saturation  modes,  the  majority  hole 
current  may  even  exceed  the  minoricy  eleccron  current.  This  implies 
chat  che  currenc  gain,  p,  could  be  smaller  chan  unity. 


From  <3, 2)  and  the  inclusion  of  position-dependent  energy-gap 
narrowing, 

En  " + SWW/W  - (l/q)34E  /8x.  (3.35) 

Substituting  (3.35)  into  (3.1)  yields  Che  expression  for  the  minority 
electron  current  in  high  injection  (n  - p) 

JN  " ^n  ,uv  >JP  + fpq»n8n/Sx  (3.36) 

where  fp  is  defined  in  (3.34).  Substituting  (3.36)  into  (3.3)  gives 

an/at  - (i/q)(/*n/pp)ajp/ax  + fpDna2n/ax2  - n/rH.  <3.37) 

where  rH  is  the  recombination  lifetime  in  high  injection.  For 
simplicity  and  conciseness,  we  shall  neglect  recombination  in  Che 
following  derivation.  Solving  (3.4)  and  (3.37)  in  Che  base  layer 


n(x,s)  - (sinh[*H(W-x)]in(0,s)  - sinh(A  x)dn(W,s))/sinh(A  W)  - p(x,s) 

(3.38) 

where  Ay  - 2s/f  Da  and  is  the  ambipolar  diffuslvity. 

Substituting  (3.38)  into  the  Laplace  Transformation  of  (3.4)  and 
integrating  the  resulting  equation  from  x - 0 to  x - W gives 
jp(o.s)  - qs( (coshAyW  - l)/AHsinhAHW]4n(0,s)  (3.39) 

in  the  forward-active  operation  and 

jp(W.s)  - qs [ (coshAyW  - l)/AHsinhAHH14n(H,s)  (3.40) 

in  the  reverse-active  operation.  Substituting  (3.38)  - (3.40)  into 


(3.36)  results 


a2i(5) 


(3.41) 


All"  AqfpDnVsch(AHW) <cosh(»HW>  ' IDn/<Dn  + D )l[cosh(»HM)  - 1J)  - A22 
A12  ” ' (A1fpDn)l||csch(AHH)  “ A21 


Up  to  n 


r discussion  has  focused  on  Including  significant 
(jp)  for  high  injection,  and  we  have  Indicated  the 
need  to  Involve  models  for  the  quasi-neutral  emitter  and  the  emitter- 
base  Junction  layers.  Even  for  hlgh-0  transistors  In  which  low 
injection  prevails  In  the  base  layer,  however,  the  forward-active 
response  to  extremely  fast  transient  excitation  may  provoke  the 
condition,  Jp  . jN,  during  a short  time  Interval.  This  will  occur  only 
near  Che  emitter  edge  <x  - 0)  of  the  space  charge  layer  where  it  feeds 
the  displacement  current  that  charges  the  emitter  Junction-layer 
capacitance  (and  changes  the  corresponding  barrier  height) . To  include 
this  possibility,  we  follow  the  approach  just  described.  Equation 
(3.35)  remains  valid  for  low  injection,  where  AEg  now  reflects  energy 
gap  changes  that  result  from  space  variations  in  chemical  composition. 
But  (3.36)  is  modified  because  in  Eq.  (3.6),  on  which  it  is  based,  the 
term,  (P/n)(pp/pn)  [jjj(x,  t)/jp(x,  C)  J , may  not  be  much  greater  than  unity 
even  though  p/n  • 1 in  low  injection.  This  violation  will  make  (3.6) 
nonlinear.  If  this  violation  occurs,  however,  it  will  occur  only  near 
x - 0 and  only  for  a short  time  interval.  Inclusion  of  the  emitter  and 
base  resistances  and  models  for  the  emitter-base  junction-layer 

(as  well  as  the  model  for  the  emitter  layer),  however. 
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will  define  jp(0,t)  in  che  shore  time  interval  of  interest  when  j p > 
jfl,  thus  avoiding  che  nonlinearity.  The  analysis  then  proceeds 
straightforwardly  to  che  low- injection  replacements  for  (3.37)  - (3.40) 
and  che  corresponding  base  layer  model. 


In  Fig.  3.5(a),  we  plot  thi 
forward  collector  current  veri 
<VgE)  for  different  doping  rat 
volts) , the  forward  collector 
and  is  independent  of  doping  i 
volts) , the  collector  current 


large-signal  model  prediction  of  the 
is  the  emitcer/base  junction  voltage 
.os  K.  In  high  injection  (VBE  >0.8 
:urrenc  has  a slope  proportional  to  q/2kT 
:Cio  K.  In  low  injection  (VgE  < 0.5 
as  slope  proportional  to  q/kT  and 
prediction  differs  from  that  of  [Gl]  because  we  do 
ic  base  storoge  charge  is  independent  of  K in  che 


simulation.  In  Fig.  3.5(b),  we  si 
collector  current  plotted  versus 
except  that  the  slope  of  reverse 

e injection  level  is  st 


collector/base  junction  voltage  (Vgc) , 
to  q/kT  as  predicted  by  [61).  This  is 


In  Fig.  3.6,  we  compare  the  model  prediction  of  Ic  - VBE  against 
experimental  data.  The  transistors  measured  havo  heavily  doped 
collectors,  which  prevent  base  pushout  in  high  injection.  As  shown,  the 
model  prediction  agrees  with  experiment  in  low  Injection  and  moderate 
injection  but  disagrees  in  high  injection.  The  discrepancy  is  probably 


voltage  for  exponentially  doped  base  layers  with  different 
doping  ratio.  K - NA(0)/NA(Wb) . ^ Th  R 


„ /(AqQ,/WQ 


the  extrinsic  resistances  (Re , RB) 


>-dimensional  effects 


such  as  emitter  crowding.  Figure  3.7  compares  the  transient  collector 
current  predicted  by  the  Gummel-Poon  model  and  the  models  developed  for 
a nearly  seep  input  voltage  vBE.  The  inductance  version  of  the  all- 
injection  non-quasi-static  model  shows  intrinsic  carrier  propagation 
delay  (-  rd  by  Ashar's  definition)  in  the  simulation  while  the  Cummel- 
Poon  models  and  the  negative-capacitance  model  developed  here  follow 
Che  input  voltage  with  no  delay.  The  negative -capacitance  model,  which 
is  similar  to  other  recently  developed  models  |4|.[7],(12)  and  is  less 
accurate  than  the  inductance  model,  predicts  Che  collector  curront  has 
negative  overshoot.  This  overshoot  arises  from  large  charging  currents 
through  the  negative  capacitance. 


3.3.2  SMll-Stgpal  Mode)? 

In  Figs.  3.8(a)  and  3.8(b),  we  plot  the  forward  transit  time  versus 
Vre  and  the  reverse  transit  time  versus  VBC  respectively.  In  high 
injection,  both  Che  forward  and  the  reverse  transit  cime  are 
independent  of  doping  ratio,  and  they  all  merge  ac  one  half  of  Che  low- 
injection  values  for  a uniformly  doped  base.  In  low  Injection,  the 
forward  transit  cime  decreases  as  doping  ratio  increases,  while  the 
reverse  transit  time  tends  Co  increase.  Since  the  base  layer,  with  K - 
2089,  is  still  in  high  Injection  at  VBC  - 0.5  volts,  Che  effect  of 
doping  ratio  to  the  reverse  transit  cime  is  insignificant  in  Fig. 
3.8(b).  Figures  3.9(a)  and  3.9(b)  illustrate  the  model  prediction  of 
delay  time  versus  input  voltage.  As  shown,  the  dependence  of  delay 


Injection  levels  resembles 


(duie)  3] 


Fig.  3.6  Comparison  of  the  forward  collector  current  versus  VBE 
between  model  prediction  and  experimental  results. 


TIME/DELAY  TIME 


vo 1 Cage  derived  via  ch 
negative -capacitance  m 


esponse  to  a large-signal  seep  inp 
Gummel-Poon(GP)  model,  via  the 


The  model  predictions  of  the  transit  time  versus  input 
voltage  for  exponentially  doped  base  layers  with  different 


Tr/(Wb2/2D„) 


of  Che  crenslc  time  on  doping  redos  and  Injection  levels  in  Fig.  3.8. 
In  comparing  Figs.  3.8(a)  and  3.6(b)  with  Figs.  3.9(a)  and  3.9(b),  ve 

is  uniform  and  Che  entire  base  is  either  in  high  injection  or  low 
injection,  or  (b)  Che  voltage  applied  is  low  enough  that  Che  entire 

large  enough  that  the  entire  base  layer  is  in  high  injection.  Thus  the 

improvements,  we  compare  against  experiment  the  phase  of  h2i  (short- 
models  developed  here.  Wo  measured  large  transistors  (Wj,  - 14  pm,  A - 

measured  by  standard  s-parameter  techniques.  The  range  of  frequencies 

p A and  by  Vgc  - - 9.3  volts.  The  new  models  show  better  agreement  with 
the  experiment  chan  that  shown  by  Che  Gummel-Poon  model,  especially  for 

Figures  3.11(a)  and  3.11(b)  show  the  inaccuracy  in  base  transit  time 
and  delay  times  caused  by  excluding  Che  plasma- induced  energy-gap 
narrowing.  Since  che  built-in  electric  field  is  reduced  by  this  gap 


rd/(Wb2/6Dn) 


voltages  for  exponentially  doped  base  layers  with  different 
doping  ratios. 

(a)  The  forward  delay  time;  (b)  The  reverse  delay  time. 


Fig.  3 . 9- -continued 


predicted  if  this  gap  narrowing  is  ignored.  Similarly,  Figs.  3.12(a) 
and  3.12(b)  deal  with  the  phase  and  the  magnitude  of  hjy  (the  input 
impedance) . The  merging  of  the  phase  at  low  (f/fT  < 10‘3)  or  high 
(f/fT  > 1)  frequencies  occurs  because  the  base  diffusion  capacitance 
has  negligible  influence.  The  frequency- independent  deviation  of  the 

narrowing  decreases  the  base  resistance.  At  moderate  frequencies  (10-3 

because  this  narrowing  increases  Che  diffusion  capacitance.  Thus  the 
electron-hole  plasma  slows  the  transient  response  and  it  adds  to  the 


The  one-dimensional  non-quasi-i 
junction  transistors  developed  hel 
large-  and  small-signal  models 
for  the  first  time,  Che  influence 
narrowing  that  arises  at  high  collector 
interactions  in  the  base.  The  models 

generally  for  high  injection. 

The  non-quasi-static  large-signal  mode: 
delay  of  the  collector  current  larger  that 
models.  The  non-quasi-static  small-signal 


oodel  predicts  more  phase 


(quasi-static) 
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Fig. 


The  effect  of  electron-hole  plasma  to  the  base  transit  time, 
(a)  The  forward  base  transit  time;  (b)  The  reverse  base 


Fig.  3.12  The  effects  of  electron-hole  plasma  to  the  input  Impedance 
h; j . (a)  The  phase;  (b)  The  magnitude. 


narrowing  slows  Che  base  transit  and  delay  times,  and  adds  phase  shift 

models  Co  including  base  pushout  and  two-dimensional  effects  such  as 
emitter  crowding  will  be  discussed  in  next  chapters. 


CHAPTER  FOUR 

HIGH-INJECTION  NON -QUASI -STATIC  COMPACT  MODELS  FOR  BIPOLAR  JUNCTION 
TRANSISTORS  HAVING  LIGHTLY  DOPED  COLLECTORS 

A.  1 Introduction 

Modern  bipolar  junction  transistors  have  submicron  base  layers 
with  doping  concentration  up  to  1018  cm'3  to  prevent  punchthrough  and 
collector  doping  as  low  as  1038  cm*3  to  assure  reasonable  breakdown 
voltage  [69],  Scaling  down  the  transistor  geometry  increases  speed.  But 
this  also  increases  current  density,  which  may  result  in  high  Injection 
effects  such  as  Webster  effect  [19],  Kirk  effect  (base  pushout)  [24], 
energy  gap  narrowing  due  to  carrier-carrier  interaction  [20J - [23] , 
current  crowding  [16] , [26] , [27]  and  spreading  [29],  High  injection 
affects  the  performance  of  transistors  [16] . [19) , [21] - (24) , [26] , 

[27], 1 29].  The  use  of  a polysilicon  technology  enhances  the  device 
performance.  But  the  resulting  polysillcon-emitter  layer  modifies  the 
terminal  performance  [69].  Heavy  doping  effects  in  the  emitter  as  well 
as  base  layers  (such  as  Fermi-Dirac  statistics,  position-dependent 
diffusivlty,  lifetime  and  energy  gap  narrowing  due  to  dopant-carrier 
interaction)  need  to  be  considered  in  determining  the  performance  of 
transistors  [14], [69]. 

The  effects  of  heavy  doping  concentration  and  high  current  density 
on  the  performance  of  transistors  have  been  discussed  [14],  [16],  [19], 
[21] -[24], [26], [27], [29], [69],  The  non-quasi-static  effects,  which  take 
into  account  carrier  propagation  delay  in  quasi -neutral  layers,  become 


important  when  base  pushouc  occurs  or  when  transistors  have  poly- 
silicon emitter  layers  [14J.  For  compound  semiconductor  heeerojunction 
bipolar  transistors  (HBT) , the  non-quasi-static  effects  are  even  more 
important  because  the  base  resistance  and  the  collector  capacitance  are 
substantially  smaller  than  the  values  they  have  in  silicon  transistors 
[70). 

In  this  chapter,  we  shall  focus  on  integrating  high  injection  and 
heavy  doping  effects  into  a non-quasi-static  lumped  circuit  model  by 
using  a systematic  method  presented  before  (14).  Non-quasi-static 
circuit  models  for  heavily  doped  quasineutral  emitter  layers  in  low 
injection  and  for  lightly  doped  quasineutral  base  layers  in  moderate 
and  high  injection  have  been  developed  in  Chapters  two  and  three  and  in 
(14)  and  (15)  for  n+/p/n+  (or  p+/n/p+)  transistors  in  which  base 
pushout  is  negligible.  However,  a comprehensive  simulation  of  modern 
transistors  must  include  models  for  lightly  doped  collector  layer  in 
which  base  pushout  occurs  at  moderate  or  high  injection. 

In  moderately  doped  base  layers  (NA  - 1018  cm'3),  AE-  (energy  gap 
narrowing)  = 50  mV  is  predicted  at  low  injection  (71)  and  AEg  < 50  mV 
is  predicted  at  moderate  injection  (22).  Energy  gap  narrowing  at  low 
injection  results  from  dopant-carrier  interaction  only.  But,  at  high 
injection,  the  dopant -carrier  interaction  decreases  because  of  the 
screening  caused  by  excess  carriers.  Thus  energy  gap  narrowing  in  the 
base  decreases  (22). 

In  lightly  doped  collector  layers,  the  carrier  concentration  can 
be  high  enough  to  decrease  the  energy  gap  because  of  carrier- carrier 


density)  - 


interaction  or  Debye  screening.  For  An  (excess  carrier 
5xl017  cm"3.  AEg  = AO  mV  [4). 

In  Section  4.2,  we  derive  non-quasi-static  models  including  energy 
gap  narrowing  for  heavily  doped  quasineutral  base  layers  at  low  and 
high  injection  by  applying  Che  method  presented  in  Chapter  two  and 
[14].  Me  Chen  generalize  the  models  Co  all  injection  levels  by  using 
Che  fitting  scheme  in  Chapter  three  and  in  [15J  and  [61).  In  Section 
4.3,  following  a short  description  of  the  mechanisms  occurring  in  the 
collector  layer,  we  extend  Che  non- quasi -static  models  to  Include  Che 
Kirk  effect,  or  base  pushout,  in  the  collector  layer.  We  show  the  model 
prediction  of  the  transit  and  delay  time  at  various  input  voltages  in 
Section  4.4  and  conclude  this  work  in  Section  4.5.  In  the  derivation, 
we  consider  a one-dimensional  n**/p+/n  transistor  shown  in  Fig.  4.1. 

4.2  Models  for  Base  I-avers 
4.2.1  governing  Equations  for  Model  Derivation 

As  shown  in  Soccion  2.2,  the  one-dimensional  carrier  transport  in 
quasi-neutral  layors  of  bipolar  junction  transistors  is  described  by 
two  current  equations  and  two  continuity  equations 
J„<x,t)  - qpn(x)n(x,t)En(x,t)  + qDn(x)On(x,t)/8x)  (4.1) 

jp(x.t)  - qPp(x)p(x, c)Ep(x, t)  - qDp(x) (3p(x, t)/3x)  (4.2) 

3n(x,c)/3t  - (l/q)«JN(x,t)/ax)  - n(x,c)/Tn(x)  (4.3) 


3p(x,t)/at  - -(l/q)(3jp(x,t)/3x)  - p(x,t)/f 


(SCL)  N (QNC) 


'] 


Wb  xb  xc  = Wb+Wc 


Charge  Layer)/QNC(Quasi-Neutral  Collector)/contact  and  P 
QNB/QNC/SCL/contact , are  shown  for  the  collector  loyer. 


toy  (l/2>(  [J^j*  *n*e(Mtp<flVnFJ/W  J-l)]1/Z  - NAJ)  j-1,2....  (4.5) 

Hero,  Vqpj (x  - xj  ) - vpj  - vgj  is  the  difference  between  the  hole  and 
electron  quasi-Fermi  potentials,  and  n^e  is  effective  intrinsic  doping 


(4.6) 


where  AE-  is  apparent  energy  gap  narrowing,  which  includes  Fermi-Dirac 
cm'3  [20]. 


4. 2. 1.1  Governing  equations  for  hleh-iniactlon  models 

injection  (Ap  = An  » NA) , the  electric  field  acting  on  holes,  Ep  - 
d(Ev/q)/dx,  differs  from  the  electric  field  acting  on  electrons,  En  - 
d(Ec/q)/dx,  because  of  energy  gap  narrowing 

En(x,t)  - Ep(x,t)  - (l/q)3AEg/3x  (4.7) 

where  Ey  and  Ec  are  valence-  and  conduction-band  edges  respectively. 

interaction  or  both  [22].  From  (4.2), 

Ep(x,t)  - (jp  + qDpSp/axj/qSpP.  (4.8) 


By  substituting  (6.7) 


and  (4.8)  into  (4.1),  wo 

jN  - qDn( (|inDp/apDn+  l)3An/3x  - <>■„/»„) (l/q)4n8AE g/3x]  + (nn/np)3p- 

(4.9) 

The  Einstein  relations  for  a heavily  doped  p-type  base  layer  in  high 

Dn/pn  - kT/q  and  Dp/Mp  - (kT/q)  F1/2<’'fP)/F-1/2(,,FP  > <*•»» 

where  qpp  - (Ey  - Epp)/UT  and  Fj/2  is  Fermi-Dirac  integral  of  order  of 
1/2.  Substituting  (4.10)  into  (4.9)  yields 

jN(x.t)  - qDnI(F  + 1)  3An/3x  - (l/kT)an  dAE  /3x]  + (f*n/u  )jp.  (4.11) 


where  F is  defined  as  Fl/2/F-l/2'  "hlch  can  be  approximated  by  a 

(71),  by  an  analytic  approach  [21]  or  by  an  numerical  approach  [22] . 
ln(An) , we  may  fit  AEg  as  a piecewise  linear  function  of  ln(An) 


= a(kT/q)ln(n/n 


» a(kT/q)ln(An/n,  ; 


4.12) 


Substituting 
- qDn(F  + 1 


by  AEg(n). 


qD  f„3An/3x 


- a/q)3An/3x 


i.13) 


Substituting  <4.12) 


explicitly  show  the  dependence  of  An  on  vqF 

A"j-  |nlo/exp(b/kT)l(exp(qvQFjAT>  - 1).  (4.14) 


4. 2. 1.2  gpvprnlnt  eqpqtlpns  for  low- inlfction  models 

For  low  injection  (An  = Ap  « HA) , by  assuming  nondegeneracy  for 

En(x,t)  « En(x)  * Eno(x)  - (kT/q)d(ln  NAof£)/dx  (4.13) 

where  F-n0  is  the  electric  field  acting  on  electrons  in  equilibrium  and 

Substituting  (4. IS)  into  (4.1)  yields 

j„(x,t)  - qDn(x)iAn(x.t)d(lnNAcff)/dx  + dAn(x,t)/dx) . (4.17) 


In  low  Injection, 
of  ln(An) 

AEg  = a- (kT/q)ln(Nfl/n1( 
From  (4.2),  (4.7)  and  ( 


a piecewise  linear  function  of  ln(NA)  instead 


I - b'.  (4.18) 


JN  - ql>n(F/M  * l)8An/8x  + (an/a  H)jp-  qDnfLSAn/ax  + (|i*«)jp  (4.19) 


W“  - C*A«i'W  I/V"'1 
*!!<-■■>  - ^ VWW  J7  «V«»  * '<*!>“!  ■ 

- J%v*1MV*iiW  J/V"!1  * "V*! 

to  .hito,  to,  bito  toi—to.,  I,  - IH,  , - v>,  ■“ 

A2(x,s)  - a + l/rh(x)  and  a - d/de; 


A2(x,s)  -st  l/rn(x) . 

Here,  rh  and  rn  are  recombination  lifetime  in  high  and  low  Injection 
respectively.  Since  the  base  recombination  lifetime  is  usually  smaller 
than  the  base  transit  time,  we  shall  neglect  the  recombination  in  the 
following  derivation.  From  oqs.  (4.4)  and  (4.20),  we  obtain 

jp(0,s)  - qAn(0,s)  I Jgb|  1 + E B21(x,s))dx  - J0bE  AjjjCx.sjdx- 

(1  + E B21(Hb,s)|/S  *21.1(»b.s)>/llA2  + J0bS  C21  ^(x.sjdx 

- fnbE  A,,  ,(x,s)dx-[E  C„  (W  s)/EA  (W  s))| 

JO  Zi  1 il  1 b Z1  1 b (4 ,21) 


If  we  integrate  (< 

A-j„(0.s>  - i1(s),  A-j, 
«b> , A"!  - An(x  - 0)  a: 


(4.3) , (4.4),  (4.13)  and  (4.19)  from  x to 
*e  obtain  jg(Ub,s)  and  jp(Ub,s).  By  defining 
0 ■ - i2(s) , v2  - vqf(x  - 0),  v2  - vqF(x  - 
; ■ An(x  - Mg) , we  obtain  the  following  two- 


1,U) 

i2(s) 


(4.22)  to  the  first  order  of 


in(s)  - [Au(s>  + A12(s)]An1(s) 


hi(s)(T) 


(j)  hi  (s) 


Fig.  4,2  A network  representation  of  the  two-port  equation  (including 
the  majority  hole  current) . 


[A22(s)  + A21(s)]An2(s) 


I-A12(s)1  [An^(S)  - An2(s)) 

1(h,X,<1  + "dft.iJ1"1  |inl(s)  ’ in2(s)1  <4'25) 

ft* * " *>* * Cv; “ ‘‘‘ ■ 

C-fe  O1  * " Cvs  '-C%  *■  '•■**> 
Cf-w^'C-4" 


(4.30) 


C11  " *11"  ’i(h.l)  /C11 

C22  " I<h.i;fTr0i.j;,,4n2{t)^v21,  R22_  ’i(h.l)  /C22 
R - (V1  - »2  “,,)].  *•  - 'd(h.i;R 


(4.25),  Co  Che  circulc  model,  (4.30) - (4. 32) . Decells  of  a similar 
derivacion  are  in  [14]  and  |15]. 


3 Hgdela_f|g. 


Based  on  Che  high-  and  low- injection  models  derived  in  SecCion 
4.2.2,  we  now  develop  non-quasl-static  models  for  all-level  injeccions 
by  utilizing  Che  fleeing  scheme  proposed  in  [IS]  and  [61].  All- 
injeccion  models  have  Che  same  circulc  configuracion  as  high-  and  low- 
injecclon  models  except  chac  Che  circulc  elemencs  contain  fleeing 
parameters.  The  circulc  elemencs  of  all-injection  model  for  large- 


c22  “ TR  IR-[4n2(C)/v2J,  Rj2-  TDR/C22  (4.34) 
R - (vt  - Vj  J/IIF-Aly  IR-4n2)  (4.35) 
L * (W1  " v2  /™F  - IR-4n2/TDR]  (4.36) 
where  Che  expressions  for  the  fitting  parameters,  IF,  IR,  TR,  TF,  TDF 


Fig.  4.3  The  large-signal  quasi-scacic  model  for  the  base  layer 
(including  the  majority  hole  current). 


For  small-signal  excicaclon,  ve  assume  exp(qvj/kT)  a exp(qVj/kT) (1 
+ qfvj/kT),  follow  Che  derivation  of  Seccion  4.2.2,  and  obcain  Che 

Cu  - TF  If-fqiNjAT),  ru-  TDF/cu  (4.37) 
c22  - TR- Ir- (qANg/kT  ),  r22-  TDR/cJ2  (4.38) 
r - kT/(qIr-4N2)  , Lss  - TDR/r  (4.39) 
6n  - q(If-AK1  - Ir-AN2)/kT  (4.40) 
where  Che  expressions  for  Che  ficcing  parameters,  If  and  Ir  are  also 


4.2.4  Remarks 

The  physical  meaning  of  Che  circuit  elements  — Ru,  R2j,  On, 

C22,  R and  L — has  been  discussed  in  [14]  and  [15].  The  majority  hole 

only  when  the  current  gain  (p)  is  small,  a condition  Chat  prevails,  for 
example,  for  transistors  operacing  in  Che  reverse-active,  quasi- 
saturation  and  saturation  modes  or  at  high  injection. 

Fitting  parameters  listed  in  Appendix  D explicitly  show  the 

neutral  emitter  and  collector  layers  into  Fig.  4.3  will  yield  a 
complete  one-dimensional  large-signal  model  as  shown  in  Fig.  4.5.  The 
expressions  for  Re,  Ce  and  Rod,  which  are  associated  with  the  injected 
hole  curronc  into  the  quasi-neutral  emitter  layer,  are  presented  in 
[14]  and  Rc,  Cc  and  Rcd , which  are  associated  with  the  injected  hole 


derived 


Fig.  4.4  The  small-signal  quasi -static  model  for  the  base  layer 
(Including  the  majority  hole  current). 


Che  approximations  for  the  Fermi-Dirac  integral,  and  on  the  empirical 
models  used  for  mobility,  lifetime  and,  especially,  energy-gap 
narrowing.  Since  the  empirical  models  for  energy-gap  narrowing  in 
(4.12)  and  (4.18)  cannot  cover  the  whole  application  range  (1017  cm-3  < 
N,  An  < 107°  cm'3),  we  use  different  models  for  energy-gap  narrowing  in 
different  ranges  to  achieve  better  accuracy.  For  example,  according  to 
data  measured  in  (71],  a'  in  (4.18)  is  approximately  0.61-q  and  b'  - 


An  > 1019  cm'3  . Similarly,  a and  b in  (4.12)  ere  determined  for 
different  An  in  (IS). 

The  low- injection  model  developed  in  Section  4.2.2  differs  from 
chat  in  | 14]  because  of  the  use  of  empirical  model  for  energy- gap 
narrowing  and  the  inclusion  of  the  majority  hole  current  in  model 
derivation. 


4.3  Models  for  Collector  Lavers 

4.3.1  Overview 

A lightly  doped  collector  layer  is  widely  used  to  assures  high 
base/collector  breakdown  voltage  and  low  base/collector  junction 
capacitance.  However,  because  the  electron  velocity  cannot  exceed  the 
scatter- limited  velocity,  electrons  will  accumulate  in  the 
base/collector  junction  layer  in  greater  numbers  as  the  collector 

distribution  of  Che  electric  field  in  the  junction  layer  and  will 
expand  this  space-charge  layer  (SCL)  [24] , [74] . [75] . If  we  keep 


The  large-signal  non-quasi-static  model  for  a n**/p+/n 

current).  CBE  and  Cgg  are  the  space-charge  capacitances  in 
the  base/eroltter  and  base/collector  junctions. 


screeching  (base  pushouC)  [24] , [74] , [78J  and  2)  lateral  currenc 
spreading  [77],  [79).  On  Che  ocher  hand,  Che  ohmic  voltage  drop  across 
Che  quasi -neutral  collector  layer  (QNL)  tends  Co  reduce  Che  reverse 
volcage  across  Che  space-charge  layer  and  which  narrows  the  space- 
chargo  layer  and  eventually  causes  quasl-sacuration  (75) , [80] , [81] . 

dominant  tendency  [63],  The  expansion  or  contraction  of  the  space- 
charge  layer  and  base  pushout,  will  affect  Che  terminal  performance  of 

[24) , [74) , [77) , decreasing  current  gain  0 [77)  and  causing  quasi- 
saturation  (82), [83].  Therefore,  comprehensive  and  accurate  models  for 
these  effects  are  essential  for  circuit  simulation. 

Analytical  models  Including  either  one  of  the  above  effects  or  all 
of  them  have  been  reported  [24) , [74]  - [84) , Based  on  this  previous  work 

categorixe  these  effects  into  two  groups:  1}  Che  expansion  and 
contraction  modes  which  correspond  to  QNB/SCL/QNC/contact  and  2)  Che 
base  pushout  mode  which  corresponds  to  QNB/QNC/SCL/contact.  We  have 
shown  previously  the  non-quasi-scacic  model  including  these  effects 

determine  the  corresponding  circuit  elements  by  using  a "regional 


determine  the  capacitance  model  for  space-charge  layer  in  the  collector 
layer  [85], 


A. 3. 2 Model  Derivation 

The  basic  concept  of  regional  approach  is  to  make  a detailed 
division  of  the  (Ic,  VCE)  plane  into  depletion  (QNB/SCL/QNC/contact) 
and  Injection  (QNB/QNL/SCL/concact)  modes  according  to  the  collector 
doping  concentration  Nc  and  width  (or  thickness)  Wc  [ 63 ) . A simplified 
(IC,  ^cg)  plane  is  shown  in  Fig.  4.6,  where  I ^ is  the  collector  current 
that  would  flow  if  all  carriers  in  the  base/collector  junction  moved  at 
the  scatter-limited  velocity  and  1 2 is  the  current  below  which  the 
transistor  is  operating  in  the  forward-active  mode.  According  to  [7SJ, 


where  A0  is  the  collector 


s the  saturation  o: 


(4.42) 


As  shown  in  Fig.  4.6,  currents  I[  and  I 2 divide  the  <Zq,  Vce) 
plane  into  six  regions.  In  Region  I (Iq  S I j S I2),  Che  large  value  of 
VCB  bends  to  expand  space-charge  layer  while,  in  Region  II  (Ic  s Ij  s 
Ij.) . the  small  value  of  VCE  tends  to  contract  it.  Thus  Regions  I and  II 
correspond  to  the  collector  layer  in  the  QNB/SCL/QNC/contect  mode.  In 
Region  III  (Ij  slcs  I2),  Che  collector  is  essentially  fully  depleted 
and  space-charge-limited  flow  occurs  [76].  In  Regions  IV  (I2  £ Ic  £ 


Il>-  v ”2  ill*  Ic) 


»1  S12S  lc). 


more  detailed  division 
transition  regions,  in 


be  in  the  QNB/QNC/SCL/contac t mode.  A 
<IC.  VCE)  plane  in  [63]  shows  that  sma: 
he  whole  oolleetor  layer  is  occupied 


:e-chorge  layer  or  by  quasi-neutral  layer,  exist.  But, 
because  these  transition  regions  aro  small  [ 63  ] , we  neglect  them  in  the 
present  derivation. 

In  Sections  4. 3. 2.1  and  4. 3. 2. 2,  we  shall  apply  the  regional 
the  capacitance  model  for  the  space-charge  layer  in  Che  collector  layer 


4. 3, 2.1  ONB/ONC/SCl./contact  mode 

the  collector  layer  in  the  reverse -active  or  saturation  modes  (VBC  a 
forward-active  mode  (QNB/QNC/SCL/contact) . 

Che  minority  hole  current,  lp2,  injected  into  Che  collector  layer  in 


■flexpiqVgg/kT)  - 


..43) 


operating  regions  In  Ic  versus 


pare  of  Fig.  4.5,  in  which 
Rc  " ',BO/|ICo  exp(<lvBC/kT>l 


s Che  collector 


collector  layer  (Fig.  4.1] 
space -charge  layer  and  tin 


iween  the  space -charge  layer  and  Che  quasi - 
W is  the  thickness  of  the  quasi-neutral 
Note  that  xj,  is  the  boundary  between  the 

e space-charge  layer  in  the  collector.  By 

Equation  (4.43)  is  for  low  injection  only.  We  nay  follow  the 

will  become  two-port  matrix  equations  as  (4.22).  The  small-signal 
version  of  (4.44)  and  (4.45)  are 

r.  - (kT/q)/[ICo  exp(qVBCAT)].  cc  - rce/rc.  cci-  rcd/cc,  ra-Rn 

(4.47) 

In  Regions  I and  II.  by  solving  the  Poisson's  equation  (85],  we 
obcain  the  thickness  X of  SCL  (75] 


yi  - tyy/i1  • (yij)])  ' 


obtained 


from  the  conventional  depletion  model.  Note  that  xg  u Wg  because  the 
base  doping  concentration  is  much  higher  than  the  collector  doping 


*n2(Wb)  = fin(xb)  - Ig/CqA^) 


(4.49) 


where  vg  is  the  electron  drift  velocity.  Ue  may  assume  vg  = vs  in 
(4.49)  because  of  che  high  electric  field  at  xg  [63),  [ 85 ) . However,  we 


separation  of  quasi-Fermi  potentials 
from  vBC  - Iq.Rq  because  Che  quasi-Fe 


VQp(Xg»,  which  is  the 
- xg.  Note  that  v2  differs 


flat  under  reverse  bias.  Substituting  the  resulting  An2  and  v2  into 
(4.31)  and  (4.32)  will  yield  the  base  model  in  this  case. 

To  determine  the  corresponding  collector  model,  we  need  to  know 
An(xc)  and  vqp(xc).  For  the  base/colleccor  junction  under  reverse  bias, 
we  may  assume  4n(xc)  - 0.  This  implies  chat  Rc  approaches  infinity  and 
Cc  approaches  zero.  Thus  the  collector  model  in  Fig.  4.5  vanishes. 

For  the  base/colleccor  junction  under  forward  bias  (VBC  > 0) , the 
quasi-Fermi  level  is  nearly  flat.  Thus  4n(xg)  and  An(xc)  as  well  as 
vQF<xg>  and  vqp(xc)  are  well  defined  and  so  are  the  corresponding  base 
and  collector  models. 


When  Iq  exceeds  some  critical 


following  discussion. 


relocate  high  electric  field  to  the  interface  between  epitaxial  layer 

layer  is  in  the  QNB/QNC/SCL/contact  mode,  in  which  xc  - ♦ Wc  and  U - 

0.  This  implies  that  the  collector  model  (Rc,  Cc  and  Rcd  in  Fig.  &.S) 
vanishes.  To  determine  the  base  model,  we  must  know  An2  and  x^. 

In  Region  III,  if  Ic  < Ij  + I2-  (Vc/V ) - I„,  the  collector  layer 
is  fully  depleted,  which  implies  X - Wc  and  W - 0.  Here,  Vc  corresponds 

at  which  the  collector  layer  is  fully  depleted  in  the  absence  of  Ic 


V - qN.W*/2«  (4.50) 

where  < is  the  permittivity.  If  Ic  a 1Q,  space-charge- limited  flow 
occurs  and,  in  a one-dimensional  analysis,  leads  to  base  pushout  [75], 
In  this  case,  xc  - Hg  + Wc  and 

*b  - Hb  + V1  ' 1{Io  ‘ I1)/(1C  ‘ V172’  " Ub  + Hcib  (4.51) 

where  Wci])  denotes  current -induced  base  thickness.  Since  Ic  > I J , we 

mode,  in  which  the  assumption,  vd  - vs  is  invalid.  Thus  wo  use  Ryder's 


'd<E>  (631 . (851 


(,  Ej  Is  Che  electric 


non-ohmlc  region  occurs  In  Ryder's  model.  For  Ie  < I3,  the  space- 

185],  An(xb)  - 0 for  Che  ohmic  conCacc.  For  > I3,  base  pushouc 
occurs  and  [63] 


*b 


e(l  + 9)(-f  - E2)/[qMc(an2/Ne  - 1)  + »(n'/N.-l)] 

(6.53) 


dn2  - N.  + (Nc/2B)(-f/E2  - 1)1-1  + (1  + Wlg/ly  l)/(-f/E2  -1)  ]1/2> , 

(4.54) 

n'  - Nc,  e - 20  and  B - (»/(l  + I)  KkT/q)  (qN/eg*) , (6.55) 

f - -l*/<<U.nHcI,Ac>.  (4.56) 


In  Regions  V and  VI, 
sacuraclon  or  saturation  n 


follows  [63] 


transiscor  also  operates  In  the  quasi- 
s but  with  Ic  larger  than  Ij.  Thus  we  may 
chac  in  Region  IV  to  obtain  xb  and  Ang 


" W1!’ 


where  Ej  is  the  electric  field  at  which  the  electron  drift  velocity 

Based  on  An2  and  xg,  we  can  determine  the  base  model  corresponding 
to  Regions  III  - VI  by  extending  Che  integrations  in  (4.26)-(4.29)  from 
[0,  Wb]  to  [0,  wbl  + [Wb.  xb] . In  Regions  III  - VI,  note  that  the 
collector  model  as  well  as  Rq  vanishes  because  W - 0 . 


:s  of  heavy  doping 


i heavy  doping  effects.  Since  the 

he  terminal  performance  of 


non-quasl-staclc  models  have  been  presented  in  [14],  we  show  only  tl 
effects  of  high  injection  here.  High  injection  mechanisms  include 
Webster  effect,  Kirk  effect  and  energy  gap  narrowing  due  to  carrier- 


carrier  interaction,  Ue  shall  1 

high  injections  such  os  the  falloff  of 
quasi-saturation  characteristics  observ 


asi-statlc  effects  o; 


ofV„.  The 


the  literature  |24] , [741 . [7SJ  and  [63] 
t illustrate  the  forward  transit  time  0 

simulated  has  NA  - 8xl015 


Nc  (collector 


doping  concentration)  - 1x10^-*  cm-^,  Wg  - 5.6pm  and  Wc  - 10  pm.  In 
narrowing  caused  by  the  carrier-carrier  interaction  in  a dense 

electron-hole  plasma  present  are  larger  than  those  without  the 
electron-hole  plasma  because  the  energy-gap  shrinkage  provokes  a higher 


layer;  this  causes  base  pushout  to  occur  i 
which  it  would  occur  without  the  electron- 

subs  tantially. 

almost  the  same  VBE  even  when  the  collector-base 


lower  VBE  chan  that  at 


ause  the  high  electric 
more  carriers  and  thus 

As  shown  in  Figs.  4.7 


:s.  However,  a high  reverse  bias  (VCB  - 20  V) 
to  increase  gradually  after  base  pushout 

and  4.8,  neglecting  the  electron-hole  plasma 
id  delay  times  and  a higher  onset  voltage  VBE 

.ers  (larger  values  of  VBE)  are  required  to 
• layer  at  a high  value  of  VEB  than  at  a 
slight  increases  of  VBE  con  produce  enough 
o neutralize  the  entire  collector  layer.  Thus 


aiding  electric  field  in  high  injection  (Uebster  effect),  which  also 
reduces  both  time  constants  before  base  pushout  occurs. 

transit  and  delay  times  while  the  aiding  electric  field  in  high 


Fig.  4.7  Modol  prediction  of  Che 
Vgg  at  different  Vgg's. 


forward  transit 


Input  voltage 


FORWARD  DELAY  TIME  (sec) 


Fig.  4.1 


input  voltage 


VBE  foe  base  pushouc. 


A. 3 Conclusions 


a Webster  effect,  Kirk  effect  end 
arrier-carrler  interaction; 
Ferml-Dlrac  statistics,  posicion- 


1)  high  injection  effects  s 

energy  gap  narrowing  due  to  carri 

dopant-carrier  interaction; 

velocities; 

unough  to  be  implemented  in  conventional  circuit  simulators.  For 
example,  we  may  implement  the  models  developed  into  SFICE2  by  utilizing 

to  interact  with  SPICE,  or  by  modifying  Che  source  code  of  SPICE. 

Since  Che  models  are  compatible  with  the  Gummel-Poon  model  built  into 


rward  voltage  (471, 


CHAPTER  FIVE 

NON-QUASI-STATIC  MODELS  INCLUDING  ALL  INJECTION  LEVELS  AND  DC, 
AC  AND  TRANSIENT  EMITTER  CROWDING  IN  TRANSISTORS 


For  bipolar  junction  transistors  operating  at  high  current  levels 
[25]-[37],  The  related,  so-called  emitter  crowding,  falls  into  two 


DC  emitter  crowding  results  from  transverse  majority  current  in 
the  intrinsic  base  layer.  At  high  current  levels,  a sizable  ohmic 
voltage  drop  occurs  in  Che  transverse  direction  and  leads  Co  a 
nonuniform  minority  current  density  distribution  in  Che  longitudinal 
direction.  As  a result,  dc  terminal  currents  (base,  collector  and 
emitter  currents)  will  be  less  than  is  predicted  by  a one-dimensional 

During  transients  or  for  small-signal  excitation,  the  majority 
quasineucral  base  layer  through  the  resistance  of  the  intrinsic  base 


che  dielectric  relaxation  time;  (b)  the  time  constan 
the  distributed  resistance-capacitance  charging  of  tl 


1 longitudinal  directions.  As  will  become  evident.  Che 
itant  dominates  in  governing  transient  and  ac  emitter 
transients,  this  time  constanc  affects  the  delay  seen  at 
ac  excitation,  it  affects  the  phase  shift  and 
impedances. 

15]  first  incorporated  the  effect  of  ac  emitter 
expression  for  the  base  impedance,  but  neglected  the 
:cer  crowding.  Therefore  his  model  only  applies  for 
:as  currents  so  low  Chat  dc  crowding  is  negligible.  Hauser  [26]  solved 
pertinent  ordinary  differential  equation  in  the  transverse  direction 

:tended  Houser's  model  to  include  ac  crowding.  He  also  pointed  out 


crowding  and  cannot  include  Che  effects 
derived  analytical  solutions  for  dc  and 
nonuniform  linear  distributed  model  for 


127] 


i distributed  model 


develop  a small-signal  model  for  i< 
Tang  [30]  examined  the  phenomenon  c 
e layer  by  using  t 


F1ELDAY  [31]. 


crowding  (29],  [30]  have  been  developed.  Based  on  Hauser's  work,  Lary  ec 
(33], (36), [46], 

All  previous  analytical  treatments  neglected  that  the  current 
distribution  in  the  longitudinal  direction  is  nonuniform.  They  also 

lumped  circuit  model.  These  shortcomings  are  overcome  here. 

either  the  use  of  a distributed  model  or  a device  simulator  such  as 
PISCES  or  FIELDAY,  an  analytical  approximation  is  useful  for 
determining  design  principles  and  providing  physical  insight. 

models,  which  admits  use  of  conventional  network  theory  in  the  modeling 
process  [88]. 

In  Section  5.2,  we  derive  analytical  models  for  dc,  ac  and 

[14]  to  model  crowding  for  low  injection  levels.  In  Section  5.2.3,  we 

systematic  method  to  model  high  injection.  Based  on  analytical 
solutions  for  low  and  high  injection,  we  then  determine  onset 

develop  the  large-  and  small-signal  models  for  all  injection  levels  by 


JM(*.y.O  - + J,,y<*.y.t)  ay  (5.1) 

Jp(x,y,t)  - Jpx(x.y.c)‘*x  + jpyOt.y.O  *,  (5.2) 


JNy  - 1%S„y  + qOn«n/8y 
ipj,  - WpEpy  - qDp9p/8y. 

Replacing  y in  (5.3)  and  (5.4)  wich  x w3 


vo  Che  expressions  for  j fjx 
e continuity  equations  are 


tl+/p/n  vertical  transistor  as  shown  in  Fig.  5.1.  In  the  figure,  x 
denotes  the  longitudinal  direction  and  y represents  Che  transverse 
direction.  Because  eqs.  <5 . 1> - (5 . 6 J are  nonlinear  and  coupled, 

s the  Gummel  method  and  the  Newton  method  are 
o obtain  general  solutions  (40],  In  order  to  get  closed- fo 
e essential  to  deriving  analytical  models,  we  make 
owing  assumptions  based  on  physical  and  practical 


eparacions  of  the  electron  and  hole  quasi  Fermi  levels  a 
and  (Wg,0)  are  equal  to  the  terminal  voltages  vBE  and  V; 

in  low  injection  for  a forward-ai 


The  structure  of  sn  n/p/n  vertical  bipolar  junction 
transistor.  We  define  the  transverse  direction  as  y axis  and 

of  the  intrinsic  and  extrinsic  base  layers  and  the  emitter 


A.  The  hole  end  electron  nobilities  are  constant  (\  f(x,y,t)>  in  the 
base,  and  the  energy-gap  narrowing  is  negligible. 

plasma[20)-[22]  at  high  current  density  is  negligible. 

Assumptions  2 through  7 reduce  Che  complexity  of  the  models  but  are  no! 
essential,  and  their  removal  is  discussed  below. 

assumption  at  different  positions  x for  various  bias  conditions.  Fig. 
5.2  is  a selected  result  of  our  simulations  chat  confirm  assumption  3. 


(5.3). 


CURRENT  (Amps/cm  *) 


Fig.  5.2 


PISCES  simulation  of  Jpy(0,y)  and  Jny(O.y)  in  y direction 


Since  ENy  - Epy  for  a lightly  doped  base  layer,  we  obtain  jpy  by 
substituting  Epy  Into  (5.4) 

jpy«  - qDp  p0  8(lnAp)/ay.  (5.7) 

Taking  the  Laplace  transform  of  (5.7)  and  integrating  the  resulting 
equation  with  respect  to  y from  0 to  y at  x - x'  gives 

Inverse  Laplace  transformation  of  (5.8)  leads  to  one  of  the  governing 
Ap(x\y,t)  - Ap(x',0,t)  exp[-  J jpy(x\y.c)dy/qDpp0).  (5.9) 

respect  to  y from  0 to  y at  x - x1 

Jpy(x'.y.b)  - Jp/x'.O.t)  - q A(s)  J Ap(x'  ,y,t)dy  (5.10) 

An  iterative  analytical  solution  of  (5.9)  and  (5.10)  can  be 


where  the  dots  Indicate  repetition  following  the 
earlier  in  the  equation.  If  jpy(x‘ .0,t)y/(qDppo)  * 1.  a condition  that 
prevails  for  the  modern  transistors,  (5.11)  becomes 
4 P(x'.y.t)  = 4p(x\0.t)  expt  -(l/qDpP0)jpy(x'  ,0,t)  y) 

- dp(x'.O.t)  (1  *<l/qDpp0)jpy(x' ,0,t)  y).  (5.12) 

Integrating  (5.6)  from  y - 0 to  y - L instead  of  y and  employing  the 
boundary  condition  Jpy(x',L,t)  - 0,  which  results  from  symmetry  of  the 

yx',0,c).qi(s)  j 4p(x-,y,t)dy.  (5.13) 


Substituting  (5.13)  into  (5.12)  gives 

4p(x\y,t)  - 6p(x\0,t)  (1  - 1/1V2  + Dpp0(x)/»4p(x\0,t)I.]  y).  ^ ^ 

If  L/2  • Dpp0(x)/1  4p(x-,0,t)L,  wo  can  reduce  (5.14)  Co 

4p(x’,y.t)  - Ap(x‘,0,t)[l  - m(x’,t)y)  (5.15) 


where  m(x',t)  - A 4p(x'  ,0,  t)L/DpPo-  Equation  (5.15)  illustrates  Che  y- 


e x-dependence  of  Ap , we  solve  Che  electron  curl 
4p(x,0.t)  - -exp(lnK-x/2M.)/(slnh  U) ■ [sinh  A (x-W.)  Ap(0,0,t) 


Substituting  (5.1 


Ap(x.y.t).  Note  that  (5.16)  is  derived  for  e lightly  doped  bei 

- qD  (-qAn  + 8An/8x)  (5.17; 


where  n — lnK/Wb.  Substituting  (5.15)  end  (5.16)  into  (5.17)  yields 

jNx(x,y,t)  - qDn(exp(qx/2  - rayJ/sinhA^W^)  ■ 

( (qsinhAn(x-W^)  - (1/2  + 3ray/2)  - ( 1 -my ) A^coshA^(x-W^ ) ) An^i- 

[-qsinhA  x-(l/2  +■  3my/2)  + (l-iny)A  cosh!  x]An.//K] 

" (5.18) 

where  Anx  - An(O.O.t)  - f(vBE)  end  An2  - An(Ub,0,t)  - f (vBC) . 

(ll)  and  Che  colleccor/base  junction  (12)  respectively.  The  expressions 
for  these  two  currents  can  be  put  in  the  form  of  che  following  two-porc 


(5.19) 


- (AqDn)[Ancoch(lnW|>)exp<-m^L)  +■  ( lnK/2Wb>cxp(  - 3m^L/2 ) ] 

- (AqDn)[»ncoth(AnWb)exp(-n2L)  - ( lnK/2«b)oxp ( - 3m2L/2) ] 

- - (AqDn)Ancsch(AnWb)cxp(-n1L)/yiC 

- - (AqD  )A  csch(A  Wb)exp(-m2L) -Tic 


Truncation  of  (5.19)  for  a uniformly  doped  base  with  negligible 


AqDn  1 + slr.tAnj/N^)  + r^] 

*11  V ' 1 + s((7/3)rcy(An1/lIA1)  + r^) 

1 + a[rdv(An2/NA?)  + r^) 

22-  W,  ' U.[(W)Ity(Vj!)*g 

_AqD„  1 + srdv(An1/KA1) 

"12  " Wb  ' 1 + *l<7/3>Vanl/NAl>  + 'dxl 

_A^DS  1 + Sr<jv(flng/Mft2) 

*21  ” Mg  ' 1 + =[(7/3)rty(An2/HA2)  + r^J 


(5.21) 

(5.22) 

(5.23) 


Because  An  « NA  in  low  injection,  the  contributions  of  rgy  and  rCy  in 
(5.20)-(5.23)  are  negligible.  This  Implies  that  emitter  crowding  is 


5.2.3  HtKtl  Injection 

Although  emitter  crowding  is  insignificant  for  low  injection,  ii 
can  occur  prominently  if  high  injection  prevails.  Emitter  crowding 

current,  exceeds  Ey  in  low  injection.  This  assures  validity  of  Che 

From  jpy  » jNy  or  Jny  « 0.  we  obtain  Ey  and 
jpy  - -2qD  dp/dy.  ( 


Jpx.  ' 


r Ex  and  obtain 


In  deriving  (5.26)  and  (5.25),  we  have  assumed  no  plasma* induced  energy 
gap  narrowing  present  [20 J - [ 22J . Thus  the  electric  field  acting  on 
holes  equals  that  acting  on  electrons.  An  extended  high-injection 
model  that  includes  plasma • induced  energy-gap  narrowing  has  been 


[IS]. 


presented  in  [15],  2D„  in  (5.24) 


and  2D«  in  (5.25)  should  be  replaced  with  1.3  Dn  and  1.3  Dp  if  a high 
density  electron-hole  plasma  is  present. 

Taking  the  Laplace  transform  of  (5.24)  and  (5.25)  and 
substituting  the  resulting  equations  into  Laplace  transformations  of 

variables  j px  and  dn  - dp.  Eliminating  j px  from  these  two  equations 
results  in  the  following  governing  equation  for  high  injoction 


(5.26) 


where  A^  - (1  + srb)/(DArb),  °A  ls  che  ambip°lar  dlffusivity  and  r.  is 

the  recombination  lifetime  in  high  injection. 

Application  of  the  product  method  to  (5.26)  leads  to 
dn(x.y.s)  - P(x,s)-Q(y,s).  (5.27) 

Substituting  (5.27)  into  (5.26)  yields 


uj;  - w)-p  - a2p/dx2  - Af^p  (5.28) 

(»£  + WQ  - a2Q/3y2  - A^y-Q  (5.29) 


independent  of  s,  an  assumption  that  one  can  show  introduces  negligible 
error  in  the  solutions  that  follow.  This  implies  that  v depends  on  do 
emitter  crowding  only.  Using  the  boundary  condition  jpy(x,L,s)  - 0,  we 


cosh[A.  (L-y)J 

ssh(AhyL)slnh(AhxWb)  ,tinh'A  hSV10  >A"l  + slnh<Ahx*)d"2 > • 


Substituting  (5.30) 


(5.25)  yields 


eoshlxa-y)] 

Jnx  - 2<">n  c„h(AhyL)sinhah;i») 


Substituting  (5.30)  and  (5.24)  : 
jpx  ’ 3px<0'y'a)  ' ‘*<s  * VO  I 


:"hl*hv(L-y» 


' h'  Vhxcosh<ihyL,Slnh<ihxUb) 
l[eosh[Ahlt(Wb-x)]-cosh(AhxWb)]4n1  + (coshU^xJ-llanj) . (5.32 

To  determine  the  dc  value  of  v in  the  forward-active  mode,  we 
neglect  net  recombination  in  the  emitter  junction  space-charge  layer, 
which  implies  chat  jpx(0,L)  - jg(-x60r'W , where  Jg(-xscr,L)  is  the 
minority  hole  current  density  injected  into  the  emitter  layer.  Thus 


JZco h/JjL  - 


where  JgQ  is  Che  saturation  current  density  in  the  emitter  layer  and 
Jv fWb  « 1 and  Jv fL  ■ 1,  inequalities  that  are  valid  in  low  and  moderate 


order  in  fiTt Wb  and 
to  values  obtained 


■/vfb.  Fig. 
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Che  particular  transistor  studied.  Similarly,  we  can  i 
assumption,  jpx(wb>L)  “ Jctwb+xscr*L)  • co  determine  cl 


To  obtain  the  two-port  equation  for  high  injection,  we  apply  the 
integration  scheme  in  Section  5.2.2  to  (5.31).  We  ignore  (yn/(ip) j px 
in(5.31),  which  introduces  negligible  error  L£  fi  > 10.  After 
truncation,  Che  matrix  elements  of  Che  resulting  two-port  equation  (set 

2AqDn  (1  + vfL2/6)(l  - v£w£/2) 


b (1  + w£L  /2)  (1  - vfU‘/6) 

* .[rj/d  ♦ w^/6)  + r^/q  ■ vfw^6)) 
+ >1^/(1  * -fL2/2)  e r^/U  - vfW2/6)) 

2AqDn  (1  + wrL2/6) 

Hb  (1  + w L2/2)(l  - w W2/6) 


*V(1  + 


1 + slrjy/a  + wrLV2)  + r^/d  - w^/6)] 

A22  " AU{"f  — : » »r>  (5-36) 

*21  " *12("r  — > "f*  (5-37) 

where  — > v means  o , replaced  with  v . 

In  (5.34)-(5.35),  r^y-  LZ/4Dp,  L2/12Dp,  W2/4Dn,  r^-  W2/12Dn 

describe  the  reverse-active  operation,  we  can  include  j px  in  (5.31), 
which  will  change  Che  expressions  for  Che  A parameters  in  (5.34)-(5.37) 


5.2.4  Onset  of  Emitter  Crowding 

crowding.  Although  we  shell  deal  only  with  dc  emitter  crowding  in  this 
Integrating  (5.19)  and  (5.31)  from  y - 0 to  y - L at  x - Wg  will 

the  relation,  Iq  (onset)  - 1q  (low  injection)  - Ig  (high  injection), 
which,  for  example,  for  uniformly  doped  base  and  omitter  layers,  yields 


Ic(onset)  - (2qZLDnWeND)/(2L2/3  - U2/3).  (5.38) 

The  onset  electron  density  at  (0,0)  in  Fig.  5.1  is 

An^onset)  - (qDp«bn2)/IJEo(2L2/3  - W2/3)] . (5.39) 

Equations  (5.38)  and  (5.39)  are  based  on  first-order  truncations 

onset  emitter-base  voltage  for  emitter  crowding  at  (0,0)  from  (5.39)  by 
using  the  standard  relation  between  density  and  voltage. 


The  large-signal  version  of  the  two-dimensional  non-quasi- 
statle  circuit  model  for  the  quasl-neutral  base  layer.  The 
small-signal  version  has  the  same  circuit  configuration  as 
the  large-signal  version.  But  the  expressions  for  circuit 
elements  differ. 


> 2.5  All-Inlection  Models 


In  this  section,  we  develop  all-injection  models  by  utilising  the 
the  second  Cauer  form  [891  of  the  resulting  s network  representation 


(1  + vfL2/2)(l  - vfwj;/6) 


- <v1/an1)/IP 


(■'fWb/2)(l  + -fL‘/6) 

L1-  a/Rf)l[<-  v£W2/2)TD2F  + (1  + v£L2/6)TTXF] (1  + vfL2/2)(l  - v£W2/6) 
+ (1  + vfL2/6)(vfW^/2)[(l  - wfW2/6)TTYF  + (1  + vfL2/2)TDXF] )/ 
Id  + v£L2/2)2(1  - vfW2/6)2]  (5.41) 

.1  " U/CuOVllmF/U  + vfL2/2)  + TDXF/(1  - vfW2/6)J  (5.42) 

:£  - (1  + vfL2/6)/[(l  + i/£L2/2)(1  - w£U2/6)1  (5.43) 

>S£  - ITDYF  - (1  + v£L2/6)  [TTYF/(1  + w£L2/2)  + TDXF/(1  - v£U2/6)|V 
([(-  w£U2/2)TDYF  + (v£L2/6  + 1)TTXF]  ♦ (w£W2/2)(l  + v£L2/6) 
[TTYF/(1  + w£L2/2)  ♦ TDXF/(1  - v£W2/6)  ] ] (5.44) 


where  — > v means  replaced  with  v , e 

In  eqs.  (5.40)-(5.49).  TDYF  and  TDYR  a 
associated  with  delay  in  y direction  fi 
reverse-active  modes.  TDXF  and  TDXR  del 


the  delay  time  In  1 


direction  for  i 


transit  tine  i 
expressions  fo 


to  arbitrarily  and  heavily  doped 
Che  fitting  scheme  that  leads  Co 


for  Che  forward-active  and 

these  fitting  parameters  (IF,  IR,  TTXF, 
1YF  and  TDYR)  are  listed  in  Appendix  E. 
fitting  parameters  for  exponentially 
the  fitting  parameters  for  arbitrarily 

(F  and  TTYR  in  Appendix  E still  apply 
ise  layers.  Detailed  descriptions  of 
te  parameters  listed  in  Appendix  E 


appear  in  [15).  161) . 


The  only  differences  are  Che 


1 of  Che  two-dimensional  circuit  model 
is  Che  large-signal  version  in  Fig.  5. A. 
expressions  for  the  circuit  elements.  If 
Lme  and  the  delay  time  for  large-signal 


1 by  replacing  v^Cvg)  , 


direction,  w 

-th  kT/q,  If  and  Ir.  The  expressions  for  the 

Chapters  Two  and  Three,  two  additional  voltage- 

■D  model  due  to  the  effects  of 

Che  physical  meaning  of  Che 


of  dc  crowding,  appear  in  the  2-D  model 
2-D  model  differ  from  Chose  in  the  1-D  r 
current  crowding.  Detailed  descriptions 
circuit  elements  in  Fig.  5.4  (except  Rj 

The  resulting  (2-D)  circuit  model  i 
crowding  for  a transistor  operating  in  t 
mode.  Since  Che  doping  profiles  in  emitt 

non-reciprocal.  Therefore,  in  the  2-D  a 


i [14]  a 


1 [15]. 


le  reverse-active  or  saturation 
sr  and  collector  layers  differ, 


e-active  modes  individually  and  split  the  current 
in  Fig.  5.4  into  four  current  sources  instead  of 
e 1-D  model.  They  are  1)  forward  quasi-static 


current  (vj-QSf/Rf),  2)  reverse  quasl-sceclc  current  (v2qsr/Rr),  3) 


In  this  section,  we  ess 
excitations,  namely,  dc,  ac 
effoct  of  emitter  crowding  o 


with  experiments.  The 


3 Assessments 

d transient.  Me  first  demonstrate  the 

r.  This  demonstration  is  based  on  the 
in  the  model  derivation.  We  then  show  thi 
dels  and  compare  the  model  predictions 


s.3.1  pc  crowding 

Based  on  the  analytical  solutions  for  the  forward-active  mode,  we 
plot  Che  longitudinal  minority  electron  current  density,  JuxiO.y), 
versus  y/L  in  Fig.  5.5  and  the  transverse  majority  hole  current 
density,  Jpy(x,0),  versus  x/Wg  in  Fig.  5.6.  These  figures  show  chat 
these  current  densitios  crowd  coward  the  emitter  edges  as  the  emitter- 
base  voltage  increases.  As  shown  in  both  figures,  current  crowding 
becomes  significant  ac  high  injection  levels  (Vgg  > 0.85  volt).  Fig. 
5.5  shows  that  Ji(x<0.y>  peaks  at  the  emitter  edge  (y  - 0) . Fig.  5.6 
shows  that  Jpy(x,0)  peaks  at  Che  emicter/base  junction  (x  - 0).  Fig. 
5.7  through  5.9  show  JpxiO.yJ/dgo'  n(0,y)  and  n(x,0}.  JEo  is  the 


(giuo/sduiv)  (A‘o)  XNr 


*io3 


intrinsic  base  layer  at  VBE  - 0.8,  0.85,  0.9  and  0.95  volt. 


Jpv  (x,0)  (Amps/cm2) 


xio3 


Analytical  solutions  for  the  transverse  component  of  the 
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XlO15 


Analytical  solutions  for  the  longitudinal  component  of  the 
normalized  majority  hole  current  density  in  the  transverse 
direction  of  the  intrinsic  base  layer  at  Vgg  - 0.8,  0.85, 
0.9  and  0.95  volt.  Jg0  is  the  saturation  current  in  the 
emitter/base  junction. 


n(0,y)  (l/cm3) 


xio,B 


Analytical  solutions  for  the  distribution  of  the  excess 


n(x,0)  (1/cm3) 


xio18 


Fig.  5.9  Analytical  solutiona  for  the  distribution  of  the  excess 
electron  density  in  the  longitudinal  direction  of  the 
intrinsic  base  layer  at  VBE  -0.8,  0.85,  0.9  and  0.95  volt. 


Fig. 


with  L/Ub  - 1.75  and  17.5.  In  comparison,  we  assume  chat  the 


saturation  current  in  the  emitter/base  junction.  In  Figs.  5.8  and  5.9, 

fields,  are  determined  by  the  sign  of  i/f  in  eqs.  (5.28)  and  (5.29). 

To  see  how  dc  emitter  crowding  affects  terminal  currents,  we  plot 


different  emitter  lengths  (L)  and 


o thickness  (Wb)  ratios.  As  show) 


significant  for  these  transistors  with  L/Wb  - 1.75.  For  L/Wb  - 

dramatically  from  that  of  the  one-dimensional  model  in  the  high 
injection  range  (Vgg  > 0.85  volt)  because  significant  dc  emitter 
crowding  occurs.  As  illustrated,  the  collector  current  predicted  by  tl 
two-dimensional  models  is  smaller  than  that  predicted  by  the  one- 
dimensional  model  for  L/Wb  - 17.5.  Since  the  deviation  between  the 

extract  the  value  of  intrinsic  base  resistance  from  Fig.  5.10. 

h (27],  Severson  (28]  and  Schutc-Alne  (32],  a 


lumped  representation  of  Che  base  resistance  will  inaccurately  model 
the  high  frequency  or  fast  transient  response  because  of  the 
distributed  nature  of  the  base  resistance.  Instead  of  a lumped  base 
resistance,  we  implicitly  include  Che  effects  of  the  intrinsic  base 
resistance  in  Che  circuit  elements  between  nodes  B and  E in  Fig.  5.4. 
We  use,  however,  a lumped  circuit  element  for  the  extrinsic  base 


x/Wb 


Analytical  solutions  for  the  eloctron  density  in  the 
longitudinal  direction  of  the  intrinsic  base  layer  at  VBE  - 


Fig.  5.12  Analytical  solutions  for  the  transverse  component  of  the 
hole  current  density  in  the  longitudinal  direction  of  the 
intrinsic  base  layer  at  VBE  - O.B  and  0.95  volt  and  at 
frequency  - 100  Hz  and  100  MHz. 


Che  ohmic  drop  between  c 
voltage  at  (0,0),  which 


voltage  applied  c 
e voltage  acting  o 


che  intrinsic  b, 


x/Wb  in  Fig.  5.11,  JpyCx.O)  versus  x/V 
versus  y/L  in  Fig.  5.13  for  different 

0.8  volt).  Similar  results  i 


n Fig.  5.12  and  jPlt<0.y)/JEo 
tses  and  different  frequencies. 

n(0,y)  versus  y/L  and  jNX(°.y> 

L3  illustrate  C1 


.e  strong  effect  of 

d [25)  neglected.  Since  the  current  density 
of  jpy(O.O)  at  f - 100  MHz  in  Fig.  5.12  is  due  Co  Che  crowding  of  the 


Ipx  (°.y)  /J 


y/L 


Fig. 
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14  The  one-dlmenslonal  and  two-dimensional  model  predictions  of 


Fig.  5.1 
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FREQUENCY  (Hz) 


tparlson  of  the  phase  of  hjj  predicted  by  ti 


Fig.  5.1 


disagree  more  ac  high  bias  levels  (V  - 0.95  vole)  than  at  low  bias 
levels  (V  - 0.8  volt)  because  significant  dc  current  crowding  occurs  In 
high  Injection.  The  magnitude  of  Che  input  Impedance  h^  shows  similar 


To  evaluate  the  validity  ol 


e two-dimensional  model  for  small- 


measured  has  a heavily  doped  collector  layer  which  prevents  base 
pushout  In  high  Injection.  As  illustrated,  Che  experiment  agrees  wlcl 

transistor  measured  has  L/Vf,  - 1.75,  the  Improvement  of  the  two- 
dimensional  model  Is  Insignificant  in  Fig.  5. IS,  but  Is  expected  to  1 
sizable  for  transistors  with  a higher  L/W^  ratio.  The  discrepancy 
between  the  experiment  and  the  two-dimensional  model  prediction  In  F: 
5.15  occurs  probably  because  we  neglected  base  width  modulation  and 


5.3.3  Transient  Grading 

the  Injected  minority  electron  currents  under  the  emitter  layer,  in  the 


I»x 
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i(1016 


electron  current  density  In  transverse  direction  of 
Intrinsic  base  layer  predicted  by  the  Gummel-Poon  (C 
end  by  the  two-dimensional  non-quasi -static  model  el 


transient  current  crowding.  In  the  very  beginning  <t  - 1 ps), 

spreads  through  the  whole  intrinsic  base  layer  and  reaches  steady  state 
(t  - 10  ns),  which  reflects  dc  crowding.  For  comparison,  we  also  plot 
in  Fig.  5.16  the  prediction  of  the  Gummel-Poon  model,  which  shows  no 
transient  current  crowding.  The  transient  current  crowding  illustrated 
in  Fig.  5.16  is  consistent  with  the  prediction  of  the  FIELDAY  [31J. 

5.4  Modifications  and  Extensions 


The  models  derived  in  Section  5.2  apply  to  an  intrinsic  quasi  - 
neutral  base  layer  with  a rectangular  boundary  as  shown  in  Fig.  5.1. 
When  significant  base  width  modulation  or  base  pushout  occurs,  however, 
the  boundary  of  the  intrinsic  quasi-neutral  layer  is  not  rectangular. 

For  instance,  che  voltage  drop  along  the  junctions  caused  by 
current  crowding  in  the  transverse  direction  will  result  in  an  x- 

Second,  use  Eq.  (5.30)  or  its  low  injection  counterpart,  Eq.  (5.15),  to 
determine  4n  at  discrete  values  of  y along  this  edge.  Third,  use  che 

at  each  of  these  discrete  values  of  y.  Fourth,  use  an  all-injection 
space-charge-layer  capacitance  model  to  determine  x corresponding  to 
each  discrete  value  of  y.  This  specifies  the  first-iteration  functional 


specifies 


repeated  until  the  sequence  (x(y)l  converges. 

dimensional  extended  base  thickness  when  base  pushout  occurs.  Instead 
of  using  the  separation  of  the  quasi-ferml  potentials,  we  use  the 
longitudinal  minority  current  density  along  x - Wg,  which  is  obtained 
from  Eqs.  (5.18)  and  (5.31).  as  the  initial  guess  in  iteration.  Based 
on  this  initial  guess  of  current  and  on  analytic  models  for  base 
pushout  [82) , [86] , [90] , we  can  determine  Wg(^)(y),  the 
first-iteration  base  thickness.  From  Hb<1>(y>  we  obtain  Jj&(y>.  An 

converges  to  Wg(y) , which  is  Che  base  chickness  in  the  presence  of 


To  derive  a complete 
current  spreading,  which  occ 
collector  layer,  and  the  sic 


,rs  in  the  extrinsic  be 
wall  injection  betweei 


e layer.  In  principle,  we  can  sti 
5.2  Co  derive  the  model  for  the  c 

base  layer.  Instead,  we  assume  chat  jPX  « jNX  in  tho  extrinc 

spreading.  PISCES  simulations  indicate  that  this  assumption 


injection,  we  need  to  modify  the  boundary  condition  when  we  derive  the 
current-spreading  model  for  the  extrinsic  base  layer. 


for  bipolar  junction  transistors.  The  models  developed  take  into 

emitter  crowding  and  transient  emitter  crowding.  We  developed  the 
models  based  on  the  analytical  solutions  obtained  from  solving  two- 


e intrinsic  base  layer;  end  they 

as  the  intrinsic  base 
e analytical  solutions,  we  can 
quasi-neutral  layers  in  analytic 

resulting  large-  and  small-signal  circuit  models  are  based  on  a 
truncation  scheme  which  leads  to  a hierarchy  of  Che  models  with  various 

offer  better  accuracy  than  the  Gummel-Poon  model,  we  compared  model 

The  treatment  here  assumed  silicon  bipolar  transistors.  But  the 
results  also  apply  to  hecerojunccion  bipolar  transistors  excited  by 


high  frequency  signals. 


CHAPTER  SIX 

CONCLUSION  AND  EXTENSIONS 


•e  derived  for  use  in  circuit  simulators.  The  model 
densities  and  produced  both  large*  and  small-signal  non-quasi-static 


The  models  Include  in  compact,  lumped  circuits  the  following 
aspects  of  transistor  physics: 

velocity, 

2)  position-dependent  diffusivity,  mobility  and  lifetime. 

3)  Fermi-Dirac  statistics  and  general  Einstein  relation  for 


Webster  effect,  which  is  the  aiding  effect  of  the  electric 

ulation  of  the  quasi-neutral  layers  due  to 
r expansion  of  the  space-charge  layers  under 


different 


Injection, 

9)  non-negllglble  majority  current  in  the  quasi-neutral  base 
layer  when  current  gain  (£)  is  small, 


as  well  as  in  the  saturation  and  quasi-saturation  modes. 

The  models  provide  physical  insight  and  information  for  parameter 

hierarchy,  which  enables  a tradeoff  between  model  accuracy  and 
either  by  using  user  defined  controlled  source  or  by  directly  modifying 


indicate  that  the  models  developed  have  advantages  over  the 
conventional  Gummel-Poon  model  and  other  conventional  models  in  the 

depends  strongly  on  the  accuracy  of  the  input  parameters  such  as  device 
makeup  (geometry  and  doping  profile) . The  accuracy  of  the  physical 
models  for  mobility,  diffuslviey  and  lifetime  also  affects  the  model 
accuracy. 


6.2  Extensions 

In  this  section,  we  indicate  sevetal  mechanisms,  which  we 
neglected  in  the  previous  model  derivation,  and  discuss  possible 
extensions  of  the  models  to  include  these  mechanisms. 


6.2.1  Avalanche  Breakdown 

at  high  injection,  one  can  increase  the  collector  doping  ci 

increases  and  thus  the  base/collector  breakdown  voltage  d' 

Hodels  for  avalanche  breakdown  have  been  published  [91], [92], 
them  are  based  on  a local  "lucky  electron"  model  [93] -[95],  wh 
assume  that  the  number  of  possible  ionization  per  primary  carr 
transit  is  very  large.  This  assumption  implies  the  validity  ol 
"continuous"  devices  physics  [ 96 ] - [ 98 ] . However,  in  many  modern  small 
devices,  the  number  of  possible  impact  collisions  per  primary  transit 

a finite  number  of  discrete  ionization  processed  must  be  used  to 
develop  associated  models  [99]. 


6.2.2  Collector  and  Base  Current  Spreading 

According  to  [29]  and  [79],  more  than  30»  of  the  collector  current 


spreading.  Details  of  this  extension  is  in  Section  Four  of  Chapter 


3 Current  C routline  in  the  Emitter  Laver  and  Contact 

neglecting  the  generation  and  recombination  current  inside  the 
base/emicter  junction,  we  may  extend  our  two-dimensional  models  to 


6.2.&  Models  for  Heteroiunction  Bipolar  Transistors  fHBT) 


quasi-static  models  for  heterojunction  bipolar  transistors.  Because  of 
junction.  Furthermore  velocity  overshoot  occurs  in  the  collector/base 


PROOF  OF  RECIPROCITY  IN  THE  QUASI -NEUTR 
To  prove  Ajj  - A2j  In  <2.11),  we  first  provi 
1 to  r2i  in  (2.19).  Defining  ft  - l/TVi,,  and 

■ C'>“  Jl  *W»  4“ 


,at  r12  in  (2.18)  i , 


Defining  Fj(*1)  « J f^Xj)  &r 


substituting  this  into  (A.l)  and 

the  integrals  in  both  equations  by  successively  applying  this  method, 
we  find  that  r2i  is  equal  to  T12  ■ The  high*order  terms  in  A22  or  A2 j 
con  be  shown  equal  by  Che  same  procedure. 


APPENDIX  B 

LIST  OP  THE  PARAMETERS  FOR  THE  MODELS  DERIVED  IN  CHAPTER  TWO 


- e*p(-t/r12)[(lA12)-J  exp(c'/f12)[d(CjjvJ)/dC|dt'  + 1^(0)] 

(B.l) 

- exp(-C/r12)(U/r12)-J  exp(C/r12)  lQS<C-)dC-  + 1^(0)]. 

0 (B.2) 

.(t)  -exp(-t/redM(l/rod)-J  exp(C  Aed> [dfC^J/dtJdf+l^CO)  1 . 


qAa0A0 

al/a0  - bjA0 
bl/b0 

J 6 p„(x)  [1  + B'(x)]dxAp<x>  + Spp0(Ue)|l  + BJ  <We)| 
J * P0W(1  + BJ(X)  + B(x)Ap<*)]dx  + Spp0(Mo)B(Ho) 


+ spp0<VAi(V 


1 + J 8 Pq(x)  Al(x)dx/rfi(xj 
J * p0(x)A1(x)dx 


-j; 

-r. 

-r.‘ 


dx/p0(x)Dp(x) 


tl/»p(«l)pt(s,) 
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<Ax2/2D)<2/lnK)[l  + (1  - K)/(KltlK)  | - (Ax2/2D)-f2 

- (Ax2/6D)  (6/lnK)  I (K  + 1)/(K  - 1)  - 2/lnK)  J - (Ax2/6D)f3 

- (Ax2/2D) ( - 2/lnK) I 1 + (1  - K)/lnK]  - (Ax2/2D) • E& 
uniform  doping  profile  or  for  high  injection,  f 1 . f2,  £3 


By  using  Anj/NAj 


(AqD^AxMKl+fj)  + uJ/Kl+f^fj^  + 0.5Ull) 
(AqDn/Ax)(0.5)([2/yK  + u^/L/K-fj  + 0.25u2ll 
(Bx2/2Dn)  ( [ (1+f j)  + UjJ/ld+fjj/fj  + 2u^| ) 
(Ax^D^a^/ffj-yiC)  + O.SUjl/d/fj-VK  + U21) 

- (Ax2/6Dn) I [ (1+fj)  + 0.5ul]/((l+f1)/£3  + Ul)i 

- (Ax2/6Dn)  l [ (1+f^)  + 0.5u2]/[(l+fl)/f3  + u2)| 


If  - (AqDn//lx)  ( [ (1+fj)  + UjJ/t  (1+fj)  • f1  + Ux) ) (C.  U) 

Ir  - <AqDn/Ax)(0.5)((2/VS  + Ujl/lVK-fj  + 0.5u2I).  (C.12) 

Injection  levels,  we  obtain  in(t),  i22<c)  an<*  ll2<c)  *n  c^e  lflrfie" 


l12(t> 


exp(-t/TDF)  [ (1/TDF) 
exp(-t/TDR) [ (1/TDR) 
exp(-t/TDF) [(1/TDF) 
exp(  - t/TDR)  [ (1/TDR) 


J exp(t'ADF)[d(TF-IF-An1)/dt]dt'< 
J exp(f/TDR)(d(TR-IR  An2)/dt]dt’i 
•J  exp ( t ' /TDF) IF  • An^dt  ’ + i^O)]  - 
-J  exp(t'/TDR)IR-An2dt'  + ij(0)  | . 


Ll(0)) 

(C.13) 

!2(0)1 

(C.14) 


(C.15) 


,<0)1 . 


IF  - [<!«,>  + uJ/Kl^)^  + 

tR  - (i/:fe)([ift/yic  + Uji/iVS  fj  + (1/Ih)2u2)i 

TF  - Kl+fj)  ♦ U1]/[(l+£1)/rfl+  u2Afh) 

TR  - [^/(fj-TE)  * ^1/Il/fj.yE  ♦ u2J 
TDF  - (<!+£,)  + r^l/Kl^/r^  Ujl 
TDR  - [(1+f,)  + rdhU2]/((l+£1)Ad2+  u2] 

If  - Ul+f^  + u^/Kl+f^f,  + 2Vihi 

” - tuvah/*  * * v4> 


(D.l) 

(D.2) 

(D.3) 

(D.4) 

(D.5) 

(D.6) 

(D.7) 


>.9) 


a0  - <P0/'p)  1“  + W3/(6Dprp)i  + SpcpQ(l  + W2/(2Dprp)] 
ai  - p0i“ + w3/oyp»  + Voh2/dp 
b0  - 1 + W2/(2D  t ) + S W/0 


CD. 10) 
CD.  11) 
CD. 12) 
CD. 13) 
(D.14) 
CD- 15) 


where  U is  the  thickness  of  the  quasi-neutral  collector  layer  (-  Wc  + 
Wjj  - xc)  and  Spc  is  the  surface  recombination  velocity  for  holes  at  the 
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layer  as  follows: 

IF  - (AqDn/Wb)([(l+f1)  + u1l/I(l+fl)-fl  + O.SUj])  (E.l) 

IR  - (AqDn/Wb)(0.S)([2/^  + Ujl/I^K-fj  * 0.25Uj])  <E,2) 

TTXF  - (Ub  /ED^UU+fj)  + UjJ/Kl+fjJ/fj  + 2Ul]l  (E.3) 

TTXR  - (Wb  /ED^II^/Ifj  JK)  + O.SUjJ/ll/^  M + Uj))  (E.4) 

TDXF  - <Wb  /6Dn)([(l+f1)  + O.Sl^l/Kl+fjJ/fj  + u^l  (E.5) 

TDXR  - (Wb  /6Dn)l[(l+f1)  + O.SUjjJ/Kltf^/fj  + U21)  (E.6) 

TTYF  - (L2/2Dp)  { (I/Du^U+ECI/D^/IDp/D^l) ) ) (B.7) 

TTIR  - (L2/2Dp)l(7/3)u2/[l+2(7/3)u2/(Dp/Dn+l)])  (E.8) 

TDVF-  aVeDpM^/Il+EUj/IDp/Dj+l)])  (E.9) 

TDYR  - (L2/6D  ) (u2/[1+2u2/(D  /Dn+1)])  (B.  10) 

If  - (AqDn/Wb)  ( ( (1+fj)  + u^/Kl+fj)-^  + Uj]  | (E.ll) 

Ir  - <AqDn/Wb)(0.5)((2/VE  + + O.SUj])  (E.12) 

where  fj,  f2,  f3  and  f^  are  defined  in  Appendix  C.  For  uniform  doping 


profile  or  for  high  injection,  f i , f 2 , f 3 and  f4  become  unity 
parameters  in  (E.1)-(E.12)  are  used  in  the  small-signal  model 
uj  have  to  be  replaced  by  ANj . 
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